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SUMMARY 


An  experimental  determination  of  the  Gladstone-Dale  constants 
for  dissociating  oxygen  has  been  carried  out.  The  result  obtained  for  the 
specific  refractivity  of  the  oxygen  molecule  was  in  good  agreement  with  avail¬ 
able  room- temperature  data,  indicating  that  the  Gladstone-Dale  constants 
do  not  appear  to  change  with  temperature.  The  result  obtained  for  the 
specific  refractivity  of  the  oxygen  atom  was  in  agreement  with  the  work  of 
Alpher  and  White  (Ref.  7)  but  the  present  work  was  considered  to  be  more 
reliable  in  view  of  a  much  lower  experimental  error,  and  the  removal  of  the 
assumption  of  a  value  for  the  Gladstone-Dale  constant  for  the  oxygen  molecule. 

A  theoretical  analysis  was  also  carried  out  in  which  express¬ 
ions  for  the  index  of  refraction  of  dissociating  and  ionizing  gases  were  de¬ 
rived.  These  results  were  used  in  an  analysis  of  the  formation  of  mono¬ 
chromatic  and  white-light  interference  patterns  in  a  Mach-Zehnder  interfero¬ 
meter  . 


An  appendix  is  included  giving  a  brief  summary  of  as  yet  un¬ 
completed  experimental  work  to  determine  the  Gladstone-Dale  constants  for 
ionizing  argon.  Another  appendix  gives  details  of  an  exploding  wire  light  J 

source,  which  was  developed  to  allow  the  use  of  a  Kerr  cell  unit  to  improve  f 

the  time  resolution  of  the  interferograms .  ‘ 
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1.  INTRODUCTION 


For  the  past  few  years,  one  of  the  continuing  projects  being 
undertaken  on  the  UTIAS  4  in.  x  7  in.  Hypersonic  Shock  Tube  has  been  the 
study  of  non-equilibrium  corner  expansion  flows  of  dissociating  and  ionising 
gases.  The  major  diagnostic  tool  in  this  study  has  been  optical  interfero¬ 
metry,  using  the  UTIAS  9  in.  Mach-Zehnder  interferometer,  since  such  an 
instrument  does  not  disturb  the  flow  under  observation.  However,  in  order 
to  use  this  method,  detailed  knowledge  was  required  of  the  index  of  re¬ 
fraction  of  the  gas  under  study  and  especially  the  relation  of  the  index 
of  refraction  to  the  density,  i.e.,  the  specific  refractivity,  or  the 
Gladstone-Dale  constant.  Experimental  measurements  of  these  quantities 
at  the  high  temperatures  used  in  the  expansion  studies  appear  to  have  been 
made  in  only  one  instance  (Ref.  7).  Hence  the  desirability  of  measuring 
the  specific  refractivities  to  check  the  results  of  Ref.  7.  Also,  since 
the  expansion  work  involved  a  large  range  of  $as  temperatures,  it  was  de¬ 
sired  to  investigate  the  temperature  dependence,  if  any,  of  the  Gladstone- 
Dale  constants. 


The  gases  which  are  being  used  in  the  expansion  studies  are 
dissociating  oxygen  and  ionizing  argon,  and  the  theoretical  part  of  this 
Technical  Note  is  applicable  to  both  g However,  the  experimental  work 
reported  here  concerns  the  measurement  o&.&he  Gladstone-Dale  constants  for 
dissociating  oxygen  only.  A  short  summary  of  the  work  done  to  date  on  the 
argon  case  is  given  in  Appendix  B.  This  case  will  be  treated  separately  in 
another  UTIAS  report,  when  the  work  is  completed. 


2.  THEORETICAL  CONSIDERATIONS 


In  this  section,  expressions  for  the  index  of  refraction  of 
dissociating  and  ionizing  gases  are  derived  from  Maxwell ' s  equations  and 
from  the  equations  of  motion  of  a  plasma.  A  short  discussion  of  polariza¬ 
bility  is  included  showing  how  the  index  of  refraction  of  a  neutral  gas 
arises.  The  Gladstone-Dale  constant  is  introduced  and  its  relation  to  the 
polarizability  of  the  gas  is  shown.  A  short  discussion  is  given  on  how  the 
polarizability  and  the  Gladstone-Dale  constant  can  be  expected  to  vary  with 
temperature . 

These  results  are  then  applied  to  interferometry  and  express¬ 
ions  are  derived  giving  the  intensity  distribution  of  monochromatic  and  white- 
light  interference  patterns. 

2.1  Index  of  Refraction 

2.1.1  Introduction 


It  is  well  known  that  a  light  wave  consists  of  coupled  electric 
and  magnetic  fields  which  oscillate  in  mutually  perpendicular  directions  and 
perpendicular  to  the  direction  of  propagation.  In  general,  the  electric  and 
magnetic  fields  can  be  represented  by  (Refs.  9>  10): 


if  =  *xp  [i(k*z  -  wt)] 
B*  =  B*o  exp  [i(k*z-  wt)] 


(1) 


1 


where  distance  z  is  measured  in  the  direction  of  propagation  of  the  wave, 
w  is  the  angular  frequency,  and  k*  is  the  complex  propagation  constant, 
where : 

k*  =  k  +  ifJ  (2) 

where  k  is  the  real  propagation  constant,  defined  as  27r/A  (X  is  the  wave¬ 
length)  and  P  is  the  attenuation  constant.  It  will  be  shown  (Sec.  2.1.3) 
that  S  is  related  to  the  conductivity  o-  of  the  medium  and  equals  zero  when 
a  =  0.  Thus  the  electric  and  magnetic  fields  could  also  be  given  by: 


I?  =  E?o  exp  [-  Sz]  exp  [i(kz  -  ut)]. 
B*  =  B*o  exp  [-  Sz]  exp  [i(kz  -  ut)] 


(3) 


which  shows  the  true  significance  of  the  attenuation  constant. 


The  phase  velocity,  Vp,  of  a  wave  comprised  of  only  one  fre¬ 
quency  is  defined  as  the  velocity  of  a  point  of  constant  phase  on  the  wave, 
i.e.: 

kz  -wt  =  a  constant  # 


Hence: 

or: 


* 


dz  _  _  w  s  2rv 

dt  p  ~  k  ~  277 

*X‘ 


Xv 


<*> 


where  v  is  the  frequency  of  the  wave.  In  ar  vacuum: 

Q 

Vp  =  c  *  2.9979  x  10  m/sec. 

where  c  is  the  velocity  of  light  in  a  vacuum,  a  constant  under  all  conditions. 
In  a  medium  other  than  vacuum  we  define: 


where  n  is  the  phase  index  of  refraction,  a  number  greater  than  one  in 
most  cales . . 


2.1.2  Polarizability 


(It  will  be  shown  later  (Sec.  2.1.4)  that  the  interaction  be¬ 
tween  the  medium  and  the  magnetic  field  associated  with  a  light  wave  is 
negligible.  Therefore,  in  all  that  follows,  we  consider  only  the  inter¬ 
action  between  the  medium  and  the  electric  field  associated  with  the  light 
wave.) 


If  the  medium  through  which  a  light  wave  is  passing  is  com¬ 
posed  of  electrically  neutral  "hard  sphere"  particles,  one  would  expect  that 
the  medium  would  have  no  effect  on  the  light  wave  (except  for  attenuation), 
since  there  would  then  be  no  way  in  which  the  medium  and  the  light  wave 
could  interact.  This,  however,  is  not  the  case.  An  atom,  for  example,  is 


-i*s»r 


1— - 


composed  of  a  positive  nucleus  surrounded  by  a  cloud  of  negative  electrons. 

If  there  is  no  electric  field  present,  the  centres  of  charge  of  the  nucleus 
and  the  electron  cloud  coincide.  If  an  electric  field  is  present,  the 
electron  cloud  may  be  displaced  slightly  with  respect  to  the  nucleus,  since, 
if  the  electrons  are  attracted,  the  nucleus  is  repelled,  and  vice-versa. 

This  relative  motion  of  the  electrons  and  nucleus  constitutes  an  electric 
current  which  is  present  as  long  as  the  electric  field  is  changing.  The 
steady-state  result,  when  the  electric  field  is  constant,  is  a  dipole. 

*  The  same  considerations  apply  to  molecules.  Non- polar 

molecules  behave  qualitatively  like  atoms,  as  described  above.  Polar  mole¬ 
cules  have  built-in  dipole  moments,  and  the  effect  of  the  electric  field  is 
to  rotate  the  molecule  as  well  as  to  distort  it.  This  rotation  also  pro¬ 
duces  a  component  of  relative  motion  of  the  centres  of  positive  and  negative 
charge  in  the  direction  of  the  electric  field,  again  constituting  an  electric 
current. 

In  general,  a  medium  composed  of  electrically  neutral 
particles  in  the  presence  of  an  electric  field  5?  will  develop  a  polarization 
field  P,  where  P  is  defined  as  the  induce^  dipole  moment  per  unit  volume. 

P  may  be  considered  parallel  and  directly  proportional  to  E*,  an  excellent 
approximation  for  the  electric  field  associated  with  a  light  wave,  and  for 
a  gas,  where  the  interactions  between  particles  are  completely  random,  and 
where  there  are  no  preferred  directions  suCh  as  might  exist  in  a  crystal. 

We  may  then  write: 


?=fiFisNaE  (6) 

where  p.  is  the  dipole  moment  induced  in  each  particle,  the  polarizability, 
a,  is  the  proportionality  factor  between  ?  and  E*,  and  N  is  the  number  den¬ 
sity  of  particles.  If  a  number  of  different  species  of  particles  are  pre¬ 
sent,  provided  there  are  no  inter-particle  forces  (for  example,  a  dilute 
gas  of  electrically  neutral  particles)  we  may  write: 

t  =  Z  (7) 

where  and  0!^  are  the  number  density  and  polarizability  respectively  of 
the  ith  species. 

Sc  far,  only  neutral  particles  have  been  considered.  In  the 
case  of  charged  particles,  such  as  electrons  and  ion  ,  there  is  of  course  a 
direct  coulomb  interaction  between  the  particles  and  the  electric  field 
associated  with  the  light  wave.  This  will  be  considered  separately  (Sec. 
2.1.4). 


2.1.3  The  Telegraph  Equation 

The  relations  between  the  refractive  index,  the  dielectric 
constant,  and  the  conductivity  can  be  derived  from  the  telegraph  equation, 
which  is  in  turn  derived  from  Maxwell's  equations.  The  telegraph  equation 
is  used  here  rather  than  the  more  usual  wave  equation  because  the  con¬ 
ductivity  of  the  medium  must  be  retained  in  the  analysis.  A  number  of 
quantities  must  now  be  defined  (in  rationalized  MKS  units): 


‘ 


let*  ^=€or+?, 

where  I?  is  defined  as  the  electric  displacement,  and  €  is  defined  as  the 
permittivity  of  free  space,  where 

€q  =  8.854  x  10“^  farads/meter 

If  the  medium  is  linear  and  isotropic  (i.e.,  parallel  to  and  proportional 
to  e)  we  may  write: 

if-  £0  (1  £  xe)  t 

where  xe  is  the  electric  susceptibility.  We  may  then  let: 

*e  «  (1  +  Xe) 

where  ke  is  defined  as  the  dielectric  constant.  It  is  then  noted  that: 

a  =  e  x  » 

0  e 

where  a  is  defined  as  the  polarizability,  as  mentioned  above. 

f-n0  (if  +  9), 

where  B* is  defined,  strictly,  as  the  magnetic  induction,  and  i?  is  defined 
as  the  magnetic  field,  i?  is  defined  as  the  magnetic  dipole  moment  per  unit 
volume,  or  the  magnetization.  In  this  development,  J?  will  be  set  equal  to 
zero,  since  we  are  dealing  with  a  gaseous  medium  where  the  magnetization  is 
negligible.  is  defined  as  the  magnetic  permeability  of  free  space, 

where:  „  „ 

Hq  =  kv  x  10”  henrys/meter. 

Generally,  for  non- ferromagnetic  media,  we  may  write: 

where  x  is  defined  as  the  magnetic  susceptibility.  We  may  then  write: 
m 

k*  =  <1  +  v> 

where  k_  is  defined  as  the  relative  permeability.  It  is  noted  that  the 
assumption  of  i?  =  0  is  equivalent  to  putting  km  =  1. 

It  is  noted  that  i?  is  analogous  to and  that  km  is  analo¬ 
gous  to  kg.  It  is  also  noted  that  since  k_  =  1,  B  and  if  can  be,  and  often 
are,  both  loosely  referred  to  as  the  magnetic  field,  since  they  differ  by 
only  a  constant  factor. 

We  can  now  proceed  to  derive  and  solve  the  telegraph  equa¬ 
tion.  Maxwell's  equations  are,  in  stationary  form  (Refs.  9>  and  10 ) 

V.Tf  =  q  (8) 

true 


4 


T 


1 


v  .  t  =  0 

(9) 

V  x  t  =  || 

(10) 

V  *  9  ■  ?tn.e  +  If 

. f  .  af.,  ar 

^true  <5t  €oSt 

=  j*true  +  £oke 

►  (11) 

where  qtrue  is  the  true  density  of  charge  and  dtrue  t^ie  true  conduction 
current . 

In  general,  if  the  medium  is  moving  at  velocity  v,  a  term 
must  be  added  to  the  right  side  of  Eq.  11  to  account  for  the  convection 
current,  a  term  of  the  form: 


V 

This  term  takes  this  form  because,  in  addition  to  the  actual  charge  pre¬ 
sent,  there  is  a  charge  due  to  the  spatial  variation  of  the  polarization, 
given  by: 


-qj, 

We  must  also  add  terms  to  correct  the  polarization  current  in  Eq.  11  for  the 
motion  of  the  medium.  The  total  polarization  current  may  then  be  written 
(Ref.  4): 

jf  =  (tx^l  *  ?  (vS) 

The  differential  operator  used  here  is  a  general  one,  and  can  be  applied  to 
any  vector  quantity  in  a  moving  medium.  (See  for  example.  Ref.  9  in  connection 
with  Faraday's  Law).  The  last  term  on  the  right  represents  the  current  due 
to  the  spatial  change  in  the  polarization  in  the  volume  .under  observation. 

The  middle  tern  on  the  right  represents  the  polarization  current  in  the 
volume  under  consideration  due  to  perturbations  in  the  velocity,  v*. 

.  In  general,  however,  the  spatial  derivatives  are  small 
enough  that  the  velocity  would  have  to  be  very  high  before  the  terms  added 
to  account  for  the  moving  medium  become  comparable  to  the  explicit  time 
derivatives.  In  addition,  q.  is  set  equal  to  zero,  i.e.,  there  is  no 
net  charge.  Thus  all  the  additional  terms  involving  the  velocity  may  be 
neglected. 


It  is  noted  that,  strictly  speaking,  the  polarization  field 
F  and  the  true  current  jtru_  have  components  due  to  the  interaction  of  the 
moving  medium  and  the  magnetic  field  associated  with  the  light  wave.  How- 
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ever,  these  terms  can  he  neglected  also,  as  much  smaller  than  the  correspond¬ 
ing  terms  due  to  the  electric  field,  (see  Sec.  2.1.4). 

The  telegraph  equation  may  now  be  derived.  Noting  that,  as 
defined  earlier: 

li  k  i? 
o  m 

where  1%  is  set  equal  to  one,  and  that  Ohm's  Law  states  that: 

3*  =  a  r 

true 

then,  taking  the  curl  of  Eq.  10  and  substituting  from  Eq.  11,  we  have: 

7*(7xt)  =  V  x  ||  ^ 

o7kS> 

"  -  Mo  It  (/true  *  €o  ke  S  ) 

&  v 

=  -  ^o  "  -5t  -  “o  eo  ke  ^2 
Vx  (Vxf)  =  V  (7  .1)  -  T2  E* 


Eut: 


However,  the  charge  density  q^rue  has  teen  set  equal  to  zero^  and  7.?  is 
negligible  in  relation  to  the  time  derivatives,  so  that  7.E  may  be  neglect¬ 
ed,  since: 


7j>'Vue-  7-? 

Hence,  when  reduced  to  one  dimension,  we  have: 


ar 

2  "  % 


& 


w  ov  0  0  6  atc 

which  is  the  telegraph  equation.  It  can  be  shown,  that  (Ref.  10) 


(12) 


“o  £o  =  -J 
c 


so  that  Eq*  12  becomes: 

b2^ 


-  ^  a 

b  2  0 


ar 


‘  c2  St2 


=  0 


(13) 
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Noting  that  Eq.  1  is  a  solution  of  Eq.  13,  we  can  substitute 

and  obtain; 


-  (k*)2  +  + 

€0C2 

2  V 

w  ^  =  0 
c2 

(14) 

g 

We  can  then  substitute  for  (k*)  from  Eq. 
imaginary  parts,  and  solve  the  resulting 
The  roots  are: 

2,  separate  into  real  and 
simultaneous  equations  for  3  and  k 

**--»  i  *.♦*  ? 

Jw2  2  <J2 

c2  e  c2  €q2 

(15) 

jw2  k2  +  cr2 

I  c2  *  c2  €q2 

(16) 

As  we  noted  earlier,  £  =  0  when  cr  =  0 

We  also  noted  earlier  that: 


and  that: 


so  that: 


c  k 
u 


c2k2 


w6 


Substituting  from  Eq.  l6,  we  find  that: 


Op2  ■  h  ^e  +  h 


(17) 


2.1.4  Determination  of  the  Conduction  Currents 


The  determination  of  the  conduction  currents  due  to  charged 
particles,  which  we  neglected  earlier  (Sec.  2.1,2),  may  be  done  through  the 
equations  of  motion  for  electrons  and  ions.  For  electrons,  we  may  write 
(Ref.  4): 


N  e(l?  +  v  x  5?)  +  N  e(u  x  S)  -  »  m 
e  e  e  e 


v 

c 


m 


du 

dt 


(18) 


where  Ne  is  the  number  density  of  electrons,  e  is  the  electronic  charge,  v 
is  the  velocity  of  the  gas  as  a  whole,  u  is  the  velocity  of  individual 
electrons,  mg  is  the  mass  of  the  electron,  and  vc  is  the  collision  fre¬ 
quency  for  an  electron. 


s 


i 

i 

| 


We  must  now  show  that  the  forces  due  to  the  magnetic  field 
associated  with  the  light  wave  may  be  neglected,  i.e.,  that  terms  containing 
(vxb)  and  (mx  I)  may  be  dropped.  We  have  from  Eq.  10: 


V  xtf  = 


Substituting  from  Eq.  1,  and  assuming  that  the  light  wave  is  a  plane  wave, 
so  that  all  space  derivatives  vanish  except  b/bz,  this  equation  becomes 
(performing  the  differentiation) 


(x  component) 


-  i  k*  E  =  i  w  B 

y  x 


(y  component) 


i  k*  E  =  i  u  B 
x  y 


so  that: 


v 

P 


taking  the  real  part  only.  Hence,  going  back  to  Eq.  18,  the  ratio  of  the 
forces  due  to  the  electric  and  magnetic  fields  may  be  written  (for  u  perpen¬ 
dicular  to  B): 


e£ 

euB 


e  VpB 
euB 


»  1 
u 


(19) 


Hence,  the  terms  giving  forces  due  to  the  magnetic  field  may  be  neglected. 


Therefore,  Eq.  18  may  be  rewritten: 


N  ei?  -  N  m 
e  e  e 


m 


du 


e  dt 


(20) 


If  we  now  let  the  electric  field  vary  as  described  by  Eq.  1,  and  assume  that 
the  velocity  of  the  charged  particles  is  proportional  to  the  applied  electric 
field,  Eq.  20  becomes: 


and 


elf  -  m 


eE  -  m  v  u  = 
e  c 


-  i  u  m  u 
e 


5  t  _  5 £ (vc * 1  w) _ e  /  r  sr\ 

1  u  -vc  v  2  +  u2  m  (v  2+  W2)  V  c  j 

c  ev  c  (21) 


species 

21: 


An  equation  similar  to  Eq.  20  can  also  be  written  for  each 
of  ion  present  in  the  medium,  and  the  result  is  given  by,  as  in  Eq. 


(22) 


where  v',  the  collision  frequency  for  ions,  will  generally  be  different  from 

v  .  c 
c 

8 


We  can  now  write  down  expressions  for  all  the  currents  that 
appear  on  the  right  side  of  Eq.  11.  The  conduction  current  due  to  the  motion 
of  the  electrons  is  given  by; 


J.  =  Ne  ' 

"e  e£ 


m  (v  e 
e  '  c 


77^.r-S) 


(23) 


Similarly,  the  ion  current  is  given  by: 


*^ion  “  Rion  6  uion 


=  _ ^2£l! _  (vt 

[(v')2+ a.2]  Vc 


(24) 


As  noted  in  Sec.  2.1.2,  a  change  in  polarization,  of  the  medium  gives  rise 
to  a  current.  Hence,  for  the  ions,  we  have: 

_  ion 
ion  St 


°ion  ion  St 

and  similarly  for  the  neutral  species  present: 

•3*  _  n 

Jpn-  sr 

=  a  N  |jf 
n  n  ot 


(25) 


(26) 


We  can  now  sum  Eqs.  23,  24,  25  and  26,  and  thus  form  the 
quantity  that  appears  on  the  right  side  of  Eq.  11: 


rion 


dE* 

o  St 


(27) 


However,  the  ion  current  may  be  neglected  here  as  being  much  smaller  than 
the  electron  current.  Thus,  Eq.  27  becomes: 


N  e2  v  t 


V  x  H  = 


If  we  let: 


N  e 
e 


2‘  2'  +  (°ion  Nion  +  an  Rn  ’  (  2  +  u2) 


me(vc  +  «  ) 


df? 

St 


2  ”e  e 
WP  =  eome 


where  Up  is  the  plasma  frequency,  and  substitute,  we  obtain: 


where  the  contributions  to  the  polarization  current  have  been  summed  in  one 
term.  If  we  compare  Eq,  28  with  Eq.  11  written  in  the  form: 


then  we  may  put: 


(29) 


(30) 


2.1.5  Index  of  Refraction  arid-  the -Gladstone-Pale  Constant 
Earlier  it  was  shown  that  (Sec.  2.1.3): 

''  a?" 

U)2  eo2 


np  =  h  ke  +  i  \IV 


(17) 


On  substituting  from  Eqs.  29  and  30  for  cr  and  k  ,  and  asstiming  v  «  w  , 
Eq.  17  becomes,  on  simplifying:  e 


n  2  =  (i  +  _i  £  a  u  _  \ 

P  U  11  J 


(31) 


and  then,  using  the  binomial  expansion,  we  obtain  to  a  first  but  very  good 
approximation: 

“j  <32) 


n_ 


The  number  density  may  be  written  as: 


10 


so  that  Eq.  32,  becomes 


or* 


n  -1  = 
P 


i 


a. 


V1-?  Ki»i-  i  J- 


(33) 


(34) 


Equation  34  is  the  Gladstone-Dale  equation,  and: 


(35) 


is  the  Gladstone-Dale  constant,  or  the  specific  refract ivity  of  the  ith 
species . 


The  last  term  in  Eqs.  33  and  34  may  be  written: 


1  V  Ne  e2  _  Ne  eg  y? 

2  u2  87r2  v2  m  €  87T2  c2  m  e 

e  o  e  o 

which  is  a  function  of  the  number  density  of  electrons  and  the  wavelength. 
Numerically: 

*  »  g 

I  =  4.46  x  10“l6  N  X2  =  4o.o8  % 
d 

-1 

where  Ne  is  in  particles  per  cubic  meter,  ~K  is  in  meters,  and  v  is  in  sec  . 

It  is  noted,  in  Eq.  34  that  an  increase  in  the  densities 
produces  an  increase  in  the  index  of  refraction  while  an  increase  in  the 
value  of  Up2,  i.e.,  of  the  electron  number  density,  produces  a  decrease 
in  the  index  of  refract ion- and,  if  up2  is  large  enough,  the  index  of  re¬ 
fraction  may  become  less  than  one,  i.e.,  the  velocity  of  the  light  becomes 
higher  than  in  a  vacuum.  This,  however,  does  not  imply  signalling  at 
speeds  greater  than  the  speed  of  light,,  because,  in  order  to  transmit  in¬ 
formation,  the  wave  must  be  modulated  which  immediately  introduces  other 
frequencies.  The  combined  wave  travels  according  to  the  group  index  of 
refraction,  which  will  be  introduced  later,  rather  than  the  phase  index  of 
refraction. 


It  is  seen  from  Eq.  35  that  the  Gladstone-Dale  constant  for 
a  given  species  is  a  function  only  of  the  molecular  weight  and  the  polariza¬ 
bility.  Since  the  molecular  weight  can  be  considered  constant  for  the 
species  in  question,  it  follows  that  any  changes  in  the  value  of  the 
Gladstone-Dale  constant  will  enter  as  variations  in  the  polarizability  of 
the -species,  Therefore  something  should  be  said  at  this  point  about  the 
temperature  dependence  of  the  polarizability. 

It  is  noted  from  the  discussion  of  polarization  given  in 
Sec.  2.1.2  that  the  model  proposed  for  atoms  and  non-poiar  molecules  is 
electrically  symmetrical  in  all  directions  when  no  electric  field  is  pre¬ 
sent,  i.e.,  its  electrical  "appearance"  is  independent  of  the  orientation 


taken  up  by  the  particle  in  the  course  of  its  random  thermal  motion.  Hence, 
when  an  electric  field  is  applied  to  the  medium,  the  polarization  of  the 
particle  must  be  independent  of  the  particle's  orientation  relative  to  the 
applied  electric  field.  The  polarization  process  takes  place  extremely 
quickly,  and  if  the  impressed  field  is  that  of  a  light  wave,  then  the 
polarization  varies  with  the  frequency  of  light.  Hence,  changes  in  the 
orientation  of  the  particle  will  have  no  effect,  since  these  changes  occur 
at  a  frequency  much  lower  than  that  of  the  impressed  electric  field.  There¬ 
fore,  if  the  polarization  of  the  particle  is  independent  of  the  particles 
orientation,  and  is  insensitive  to  changes  in  that  orientation,  then  the 
polarization  can  be  said  to  be  independent  of  temperature,  which  is  a  mea¬ 
sure  of  the  (relatively)  low  frequency  at  which  changes  in  orientation  take 
place  (e.g.,  10^0  collisions  per  second’,  compared  to  10^-5  sec--*-,  the 
approximate  frequency  of  light). 

The  case  of  polar  molecules,  however,  is  different.  Here,  it 
is  recalled,  the  molecules  have  built-in  dipole  moments  and  the  orientation 
of  these  dipole  moments  is  thus  a  function  of  the  thermal  motion  of  the 
particles,  when  no  electric  field  is  present.  When  an  electric  field  is 
applied,  it  tends  to  orientate  the  dipoles  in  the  direction  of  the  field, 
but  this  orientation  is  opposed  by  the  thermal  motion  of  the  particles.  It 
is  seen  that  the  degree  of  orientation  of  the  dipoles  in  a  given  electric 
field  is  a  function  of  how  strong  are  the  thermal  restoring  forces  ,  i.e., 
of  the  temperature. 

The  latter  case  of  the  polar  molecule  was  included  in  order 
to  present  a  complete  picture,  for,  in  the  present  work,  the  only  molecule 
included  is  oxygen,  O2,  which  is  non-polar,  since  both  "ends"  of  the  mole¬ 
cule  are  identical.  Therefore,  the  conclusion  based  on  the  arguments  pre¬ 
sented  here  is  that  the  polarizations,  and  thus  the  Gladstone-Dale  constants, 
of  the  species  being  studied  in  this  work  should  not  change  with  temperature . 
This  conclusion  is  however  based  on  rather  simple  classical  arguments;  it  is 
possible  that  a  quantum-mechanical  treatment  might  yield  a  different  conclu¬ 
sion,  but  this  is  beyond  the  scope  of  this  thesis.  A  detailed  discussion  of 
polarization  and  inter-molecular  forces  is  given  in  Ref.  8.  However,  much 
of  the  treatment  in  Ref.  6  concerns  liquids,  and  many  of  the  effects  dis¬ 
cussed  are  negligible  in  the  case  of  a  gas. 


2.2  Application  to  Interferometry 
2.2.1  Introduction 


Earlier,  we  denoted  the  amplitude  of  the  electric  and  magnetic 
fields  associated  with  a  light  wave  by: 

I?  =  Eq  exp  [i  (k*z  -  wt)] 

B*  =  exp  [i  (k*z  -  wt)] 

27T 

where  k*  =  k  +  ip,  k=-^*  and  w  =  2w 

If  the  medium  being  considered  is  a  gas,  then  the  conductivity, 
a,  is  small  (e.g.,  10"5  m  hos/m).  Therefore,  we  can  assume  that  there  is 
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no  attenuation  of  the  light  wave,  i.e.,  that  0  =  0  and  e  =  1  ,  at  least 
for  the  distances  involved  in  interferometry.  Therefore  we  may  write: 

E*  =  E?0  exp  [27 r  i  (f  -  vt)  ] 

A 

Up  to  now,  it  has  been  convenient  to  work  with  the  exponen¬ 
tial  form  in  expressing  the  variation  of  the  quantities  associated  with  a 
light  wave.  However  for  the  purposes  of  this  section  it  is  convenient  to 
work  with  the  amplitude  of  the  light  wave,  which  is  related  to  the  electric 
field  and  can  be  described  by  the  same  function.  It  is  also  convenient  to 
use  only  the  real  part  of  the  exponential  function  and  write: 

A  =  A  cos  2v  (vt  - 

O  A 

The  imaginary  part  of  the  exponential  could  also  be  used  for  this  purpose; 
the  only  difference  would  be  that  the  cosine  is  everywhere  replaced  by  sine. 

2.2.2  The  Mach-Zehnder  Interferometer 

Figure  7  shows  the  geometrical  and  optical  layout  of  a  Mach- 
Zehnder  interferometer.  A  more  detailed  description  of  the  UTIAS  9"  inter¬ 
ferometer  is  given  in  Ref.  3,  but  a  brief  summary  of  the  principles  involved 
is  included  here.  Light  from  the  source  is  formed  into  a  parallel  beam  of 
sufficient  area  and  directed  to  splitter  1,  where  it  is  divided  into  two 
parts  by  amplitude  division.  Fart  of  the  original  beam  is  directed  to  mirror 
1  and  from  there  it  is  passed  through  the  test  section  to  splitter  2.  The 
remaining  part  of  the  beam  passes  from  splitter  1  through  the  compensating 
chamber,  an  optical  replica  of  the  test  section,  to  mirror  2  and  then  to 
splitter  2.  At  splitter  2,  the  two  parts  of  the  original  beam  are  recombined. 
Part  of  the  reconstructed  beam  passes  through  various  lenses  and  mirrors  to 
the  photographic  plate,  ground-glass  screen,  or  other  detector.  The  remain¬ 
ing  part  passes  out  of  the  system,  and  is  not  used,  although  it  could  be  if 
convenient . 


Let  us  now  assume  that  conditions  inside  the  test  section 
and  inside  the  compensating  chamber  are  identical,  and  that  the  optica!)  nath 
lengths  of  the  test  section  and  compensating  chamber  beams  are  equal,  "his 
does  not  necessarily  mean  that  the  geometrical  path  lengths  sire  equal 
(although  they  usually  are,  in  practice),  but  it  does  mean  that  the  transit 
times  of  the  two  beams,  from  separation  at  splitter  1  to  recombination  at 
splitter  2,  are  equal.  We  also  assume  that  the  plates  of  the  interferometer 
are  exactly  parallel  and  perfectly  flat,  s 6  that  the  path  length  at  one 
point  in  the  field  of  view  is  the  same  as  at  any  other  point.  Under  r’:.  se 
conditions,  the  test-section  beam  and  the  compensating  chamber  beam,  which 
were  formed  from  the  same  beam  at  splitter  1,  will  arrive  at  splitter  2 
exactly  in  phase,  and  the  resultant  amplitude  will  be  the  same  as  that  of 
the  original  beam. 

If  a  change  now  occurs  in  the  test  section  conditions,  the 
optical  path  length  seen  by  the  test  section  beam  will  also  change.  This  is 
caused  by  a  change  in  the  velocity  of  light  in  the  test  section,  i.e.,  the 
wavelength  changes  (since  the  frequency  is  fixed),  and  these  changes  show  up 
as  a  change  in  the  index  of  refraction.  Thus,  an  optical  path  difference 
between  the  test  section  and  compensating  chamber  beams  is  introduced,  given 
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%  where  h  is  the  geometric  width  of  the  test  section,  and  subscripts  1  and  2 

!  denote  initial  and  final  conditions  in  the  test  section,  respectively.  In 

terms  of  the  indices  of  refraction,  this  becomes  (dropping  the  subscript  p 
denoting  phase  index  of  refraction) 


where  Aq  is  the  wavelength  in  vacuum. 

The  resultant  amplitude  seen  at  the  photographic  plate  is 
then  simply  the  sum  of  the  amplitudes  of  the  test  section  and  compensating 
chamber  beams,  taking  into  account  the  phase  difference  introduced.  Here, 
since  we  are  observing  in  a  fixed  position  relative  to  the  source  and  inter¬ 
ferometer,  we  can  neglect  the  z  dependence  of  the  amplitude.  It  is  noted 
that  introducing  a  path  difference  has  the  same  effect  as  adding  a  phase 
constant  to  the  test-section  beam.  Thus: 


i 


k 


*\  I 


A  =  4  A  cos  (27T  vt)  +  I  A  cos  2ir  (  vt  +  — —  -  -4f-  ) 

=  A  cos  (2tt  vt  +  ~  )  cos  7r(  -  -r~  ) 

°  ”2  x  ^1  e 

=  Aq  cos  (2tt  vt  +  JSj  -  ^  )  cos  f  (nx  -  n2) 

=  A  cos  (27 r  vt  +  —  v”  )  cos  ~  (n„  -  n  )  (36) 

O  A  ^  X 

2  1  O 

This  was  obtained  using  standard  trigonometric  relations.  No  account  has 
been  taken  of  the  light  lost  at  the  second  splitter,  since  this  alters  the 
result  only  by  a  constant  factor. 

However,  the  quantity  actually  detected  by  the  photographic 
plate  is  the  intensity  of  the  light,  which  is  given  by  the  square  of  the 
amplitude.  Furthermore,  it  is  noted  that  the  time  variation  of  the  amplitude 
is  at  frequency  v,  which  is  far  too  rapid  for  the  human  eye  or  any  other 
sensor  to  detect.  Rather,  the  sensor  sees  a  time  average  of  the  intensity. 
Therefore : 

I  =  A2  =  A2  cog2  (2tt  vt  +  y  "  IT  )  COs2  T  ^n2  "  nl^ 
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1  ■  K 

1  at  -  i  [a2  cos2  S  („2  .  ^)  j  COS2(27TVt  ♦  & 

-f  )  dt 

i 

=  2T* 

A2  cos2  (n2  -  Dl)J  {|’cos(4mt  +  ff-  ^  )it 

*{*} 

1 _ 1 

HIS 

if 

kl cos2  t  (ns  -ni>  ]  {  erf sin  (lwt  +  If*  If) 

B 

and  taking  the  limit  as  T  becomes  very  large,  compared  to  — ,  we  obtain: 

I  -  =os2  (n2  1  (37) 

0 

Thus  the  only  variation  in  intensity  which  can  be  detected  is  that  due  to  the 
change  in  conditions  in  the  test  section.  Since  we  specified  that  the  path 
length  was  the  same  throughout  the  field  of  view  initially,  we  see  that  the 
intensity  will  vary  in  the  same  way  throughout  the  field,  provided  that  n2 
is  uniform  throughout  the  test  section.  When  the  argument  of  the  cosine 
equals  an  even  multiple  of  tt/2  ,  the  intensity  is  a  maximum;  when  the 
argument  equals  an  odd  multiple  of  tt/2,  the  intensity  is  zero,  and  no  light 
is  transmitted. 


Generally,  however,  conditions  will  not  change  uniformly 
throughout  the  test- section;  rather,  n2  will  vary  from  one  part  of  the  field 
to  another  and  7 rL/AQ  (n2-n1)  will  assume  various  yalues.  When  n2  varies 
smoothly  (which  it  does,  except  through  flow  discontinuities),  7TL/A  (r^-n^) 
will  also  change  smoothly,  increasing  (or  decreasing)  in  value,  passing 
through  odd  and  even  multiples  of  7r/2.  The  result  is  -a  set  of  interference 
fringes:  the  contours  where  7rL/A0  (n2-n1)  is  an  even  multiple  of  tt/2 
show  up  as  bright  lines,  and  the  contours  where  7Tl/A0  (n2-n.)  is  an  odd 
multiple  of  tt/2  show  up  as  dark  lines,  with  a  continous  variation  in  bright¬ 
ness  between  the  bright  and  dark  lines.  The  change  in  n2  between  adjacent 
bright  (or  dark)  fringes  is  fciven  by: 


or: 


r  t(n2,  +  An)  -  nx]  =  f  (n,  -  nj  ♦  , 

n  *  n 


Therefore,  if  the  index  of  refraction  is  known  on  a  given  fringe,  the  index 
of  refraction  on  any  other  fringe  can  be  determined. 
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2.2.3  The  Parallel  Fringe  Method 

Up  to  now,  it  was  assumed,  that  the  plates  of  the  interfero¬ 
meter  were  all  parallel.  We  now  modify  the  system  by  introducing  in  one  of 
the  beams  of  the  interferometer  an  optical  path  difference  which  changes 
uniformly  across  the  field  of  view.  This  is  done  in  the  compensating  cham¬ 
ber  beam  (Ref.  3)  by  rotating  mirror  2  through  a  small  angle  €  .  The  optical 
path  difference  thus  introduced  is  (in  cycles): 


where  y  is  measured  perpendicular  both  to  the  axis  of  rotation  of  the  mirror 
and  to  the  direction  of  propagation  of  the  beam,  and  where  subscript  "a" 
indicates  the  medium  in  which  the  interferometer  is  immersed,  generally  air. 

It  is  seen  that  the  effect  of  this  alteration  is  the  same 
as  if  the  index  of  refraction  in  the  test  section  had  been  made  to  vary 
uniformly  across  the  field  of  view,  and  the  result  is  straight,  uniform 
fringes.  Mathematically,  proceeding  as  before,  the  resultant  amplitude  at 
any  point  in  the  field  is  given  by: 
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Therefore: 


I  =  A2  cos2fVt  +  #  -  +^N) 

A  A  J 


COS 


vL  (  2ey  na  \ 


Here,  as  in  the  case  of  Eq.  36,  the  time  variation  of  the  intensity  is  at 
frequency  v  ,  and  is. therefore  far  too  rapid  to  be  detected.  Therefore,  tak¬ 
ing  the  time  average,  as  before: 

1  -  if 1  dt  “  f  [Ao  cos2  T  (w  "T"2)  ]/cos2  (£svt  +  f  '  \  +  tF)  dt 

0  0  0  d  X  a.  ' 
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which,  on  simplifying  and  taking  the  limit  as  T  becomes  very  large  compared  to 
l/v  ,  becomes: 

i2  2  Ik  (  -  2eyn& 
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(39) 


This  expression  is  similar  to  that  obtained  earlier  in  Eq. 
37  for  the  infinite  fringe  method,  but  now  the  intensity  is  also  a  function 
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of  position  y.  When  conditions  in  the  test  section  are  constant,  i.e.,  * 
(n2  -  n^)  is  constant  over  the  field  of  view,  the  intensity  is  governed  by 
changes  in  the  term: 

2€yna 

L 


and  the  spacing  of  the  fringes  is  given  by: 

7tLi  2eyna  t  tt  7rL  2e(y+d)na 

T  L  A  L 
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or: 


Generally,  however,  test  section  flow  conditions  will  not  be 
uniform  and  thus  the  fringes  generally  will  not  be  straight.  Referring  to 
Fig.  8,  suppose  that  two  points  in  the  field  of  view  P(xp,yp)  and  Q(x  ,y _) 
are  on  the  same  fringe,  i.e,  that:  *  ** 
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is  a  constant.  Then: 
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Equation  4o  is  basic  to  the  use  of  the  Mach-Zehnder  inter¬ 
ferometer:  it  relates  the  y-deflecuion  of  the  frinp.es,  in  terms  of  fringe 
spacing,  to  the  change  in  index  of  refraction.  Since  and  L  are  known  in 
a  given  situation  and  the  fringe  spacing  can  be  measured,  one  can  then 
determine  the  index  of  refraction  over  the  entire  field  of  view,  provided 
that  the  index  of  refraction  is  known  at  one  point. 

If  we  now  introduce  the  Gladstone-Dale  equation,  Eq.  34, 

one  can  write: 

nP  -  \  ■  ?  Ki  (pip  -  pi«)  -  f  (“Pp  -  <“r>*  ] 

which  relates  the  change  in  index  of  refraction  to  the  changes  in  the  partial 
densities  and  in  the  electron  number  density.  If  this  was  substituted  into 
Eq.  40,  it  is  noted  that  as  many  equations  as  there  are  unknowns  would  be  re¬ 
quired  to  solve  for  all  the  mass  densities  plus  the  electron  number  density. 
If  we  assume,  however,  that  only  one  species  is  present  other  than  electrons, 
we  may  write: 
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|  therefore,  having  two  equations  of  this  form,  one  could  determine  the  mass 

t  •  density  and  electron  density  at  any  point  in  the  field  of  view,  provided 

I  that  'these  quantities  were  known  at  ore  point  in  the  field. 


The  effect  of  the  electrons  is  particularly  noted. 

If  pp  >  Pq  for  constant  ,  then  y?  >  y^.  But  if  (p,p  -  pq)  is  held  con¬ 
stant  and  *  l(wp)p  -  O^pJnJ  is  allowed  to  increase,  then  (yp  -  yQ)  will  de¬ 
crease.  So  that,  if  an  ''increase  in  dens3.ty  produces  a  positive  fringe 
shift,  an  increase  in  electron  density  will  produce  a  negative  fringe  shift. 


2.2.4  White  Light  Interference 

So  far,  only  interference  between  monochromatic  light  beams 
has  been  considered.  We  must  now  consider  interference  between  beams  of 
"white  light",  i.e.,  light  containing  a  range  of  wavelengths.  This  change 
introduces  several  complications  such  as  dispersive  effects,  which  arise 
because  the  index  of  refraction  of  a  medium  is  generally  a  function  of  wave¬ 
length. 


I 


It  is  also  impossible  to  consider  interference  between  beams 
of  white  light  in  the  same  way  as  we  considered  monochromatic  interference  - 
as  an  amplitude  phenomenon.  Bather, the  intensity  distribution  in  a  white- 
light  interference  pattern  must  be  considered  as  the  superposition  of  the 
intensity  distributions  from  an  infinite  number  of  monochromatic  interferograms 
at  wavelengths  within  the  wavelength  range  of  the  white  light.  This  arises 
from  the  fact  that  nothing  can  be  said  about  the  phase  differences  between 
the  various  wavelengths  emitted  by  the  source. 

Bor  a  monochromatic  inteferogram,  the  source  of  light  may 
be  considered  as  the  set  of  oscillators  (atoms  or  molecules)  in  the  source 
which  emit  at  the  particular  wavelength  in  question.  Each  one  of  these 
oscillators  emits  for  only  a  short  time  (but  still  several  orders  of  magni¬ 
tude  longer  than  1/  v)  so  that  the  set  of  oscillators  changes  frequently. 
However,  for  a  given  set  of  oscillators,  the  light  emitted  can  be  said  to 
have  a  certain  phase  which  is  the  resultant  of  the  phases  of  all  the  indi¬ 
vidual  oscillators;  this  phfi.se  remains  constant  for  the  given  set  of 
oscillators,  since  all  the  oscillators  within  the  set  emit  at  the  same  fre¬ 
quency.  Similarly,  if  one  considers  the  set  of  oscillators  emitting  at  a 
second  frequency,  the  emitted  light;  from  these  oscillators  can  also  be 
assigned  a  resultant  phase.  But  since  the  frequencies  of  the  light  emitted 
by  these  two  sets  of  oscillators  are  different,  the  phase  difference  between 
the  two  frequencies  is  always  changing.  Likewise  for  all  the  other  wave¬ 
lengths  emitted  by  the  source. 

However,  as  stat-ed  above,  the  set  of  oscillators  emitting 
at  a  particular  wavelength  exists  for  only  a  short  time;  when  the  set  of 
oscillators  changes,  the  resultant  phase  undergoes  unpredictable  changes 
also.  Thus  in  addition  to  the  variation  in  phase  difference  due  to  the 
fact  that  the  frequencies  are  different,  there  are  also  random  changes  in 
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phase  difference  due  to  the  changing  sets  of  oscillators,  so  that  it  is 
impossible  to  predict  the  changes  in  phase  difference  mathematically. 
Therefore,  one  cannot  consider  interference  between  beams  of  light  of 
different  wavelengths  as  an  amplitude  phenomenon,  because  to  do  so  re¬ 
quires  knowledge  of  the  phase  difference  which  is  not  available.  The  best 
that  can  be  done  is  to  consider  the  energy,  or  intensity,  distribution 
produced  by  interference  at  a  given  wavelength,  about  which  something  can 
be  said,  and  then  integrate  over  the.  frequency  spectrum  to  find  the  result¬ 
ant  intensity  distribution. 

Some  authors  have  not  considered  these  points  when  integrat¬ 
ing  over  the  frequency  spectrum  to  find  the  resultant  interference  pattern. 
Attempts  have  been  made  to  work  with  the  amplitude,  and  in  attempting  to 
define  the  phase  differences,  the  tacit  assumptions  have  been  made  that,  (l), 
the  sets  of  oscillators  which  produce  each  frequency  remain  the  same,  so  that 
no  random  changes  in  phase  occur,  and,  (2),  that  the  phase  differences 
between  all  the  frequencies  are  zero  at  time  t  =  0.  This  produces  the  not 
unexpected  result  that  the  intensity  over  the  field  of  view  of  the  interfero¬ 
meter  decreases  with  time,  and  with  distance  away  from  the  zero-order  fringe, 
which  is  not  the  case.  The  intensity  in  the  field  of  view  is  constant  in 
time,  and,  while  the  peak  intensity  does  deprease  with  distance  away  from 
the  zero-order  fringe,  the  average  intensity  is  the  same  over  the  field  of 
view. 


In  view  of  these  considerations,  a  convenient  starting  point 
from  which  to  calculate  the  intensity  distribution  in  a  white-light  inter¬ 
ference  pattern  is  the  intensity  distribution  for  a  monochromatic  interfer¬ 
ence  pattern  given  by  Eq.  39: 


-  Ao  2  VLf  26yna>| 
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(39) 


Introducing  the  Gladstone-Dale  equation  (Eq.  3*0  and  assuming  only  one 
species  present  other  than  electrons,  Eq.  39  becomes: 
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where  (<*>)2  has  been  set  equal  to  zero  to  correspond  to  the  usual  physical 
situation  of  no  electronsppresent  in  the  initial  state, v/c  has  been  sub¬ 
stituted  for  l/AQ  and  kyv*  has  been  substituted  for  u>2. 


of  the  light. 


We  must  now  determine  the  spectral  intensity  distribution 
Referring  to  Fig.  9a,  we  set: 
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V1  <  V  <  V2 
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I,  =  0  for  v  <  and  v  > 


3&C: 


The  intensity  distribution  in  the  white-light  interference  pattern  is  then 
given  by: 
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We  must  now  introduce  frequency  dependent  expressions  for  the  indices  of  re¬ 
fraction  and  the  Gladstone-Dale  constant.  This  can  be  done  through  the 
Cauchy  relation  (Ref.  14): 

2  4 

n  -  1  =  A  +  Bv  +  Cv  +  . 


K  =  Ko  +  Kgv2  +  K^v4  + 


Having  thus  accounted  for  the  frequency  dependence  of  the 
indices  of  refraction  and  the  Gladstone-Dale  constant,  Eq.  43  can  be  inte¬ 
grated,  and  the  details  are  given  in  Appendix  A.  The  result  is: 

t?  r  sin  t[X1(p£  -  gl)(A„  +  A1)  -  Xgy(l  +  a;  +  x3]  (v3-v1» 
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where:  v  = 
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and  where,  as  shown  in  Appendix  A: 
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corresponds  to  the  phase  index  of  refraction  of  the  air  in 
which  the  interferometer  is  immersed. 

(l  +  A'  +  A' )  corresponds  to  the  group  index  of  refraction  of  the  air  in 
°  which  the  interferometer  is  immersed 
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A  is  the  Gladstone-Dale  constant  of  the  gas  in  the  test 

0  section,  which  is  analogous  to  the  phase  index  of  refrac¬ 

tion,  of  the  gas,  and 

(  A  +  A1 )  is  the  Gladstone-Dale  constant  of  the  gas  in  the  test 
°  section,  which  is  analogous  to  the  group  index  of  refrac¬ 

tion  of  the  gas. 

Equation  44  shows  that  a  white-light  interference  pattern 
consists  of  a  background  intensity,  A“/4,  on  which  is  superimposed  a  sinu¬ 
soidal  variation  governed  by  frequency  V ,  which  is  in  turn  modulated  by  a 
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sine  wave  governed  by  frequency  (v  -  v  ).  It  is  also  noted  that  as  (p  -  p  ) 
varies  from  zero,  the  position  of  the  zero-order  fringe  (i.e.,  the  value  of  i 
y  for  which  the  argument  of  the  cosine  is  zero),  and  the  position  of  the 
maximum  of  the  modulating  function,  or  envelope,  (i.e.,  the  value  of  y  for 
which  the  argument  of  the  sine  is  zero)  vary  from  zero  also.  However,  the 
values  of  y  which  result  in  the  arguments  of  the  sine  and  cosine  being 
equal  to  zero  are  different,  except  when  X3  (i.e.,  the  electron  number  den¬ 
sity)  an*!  (P2  -  Pi)  equal  zero.  It  is  noted  too  that  the  contrast  of  the 
fringes,  in  addition  to  being  controlled  by  the  sine  function,  decreases 
as  distance  from  the  position  of  the  maximum  of  the  envelope  increases. 

These  effects,  and  others  described  below,  are  demonstrated  in  the  numerical 
calculations  included  in  Appendix  A. 

Of  particular  interest  is  the  contribution  of  the  electrons; 
in  the  cosine  factor  this  term  is  negative,  while  in  the  sine  factor  this 
term  is  positive.  This  indicates,  as  noted  in  connection  with  Eq.  4l,  that 
the  electrons  produce  a  negative  shift  of  individual  fringes,  but  a  positive 
shift  of  the  envelope  of  the  intensity  distribution.  It  is  also  noted  that 
the  effect  of  the  electrons  is  frequency  dependent,  since  Xo  contains  a 
factor  (V)  .  Therefore,  the  effect  of  the  electrons  will  be  less  at  shorter 

wavelengths,  toward  the  violet  end  of  the  spectrum, and  greater  at  longer 
wavelengths,  toward  the  red  end  of  the  spectrum. 

It  is  noted  too  that  Eq.  44  produces  the  correct  results  in 
the  limits  as  (vg  -  v,)-*oo  and  as  (v2  -  v^)  -*0.  As  (vg  -v^)  -»0,  the 
range  of  wavelengths  becomes  more  and  more  narrow  until,  in  the  limit,  only 
one  wavelength,  v,  remains,  and  the  white-light  interferogram  becomes  a 
monochromatic  interferogram.  In  Eq.  44,  as  (v2  -  Vi)-»  0,  the  argument  of 
the. sine  approaches  zero  for  all  values  of  (pg  -  Pq) ,  y»  and  Ne,  and  since: 

—  =  1  and  \  (l  +  cos  2x)  =  cos2x  (45) 

it  is  seen  that  Eq.  44  degenerates  to  Eq.  39 >  after  allowance  is  made  for 
the  changes  in  notation. 

As  (v2  -  v1)  -*oo  ,  it  is  seen  that  the  coefficient  of  the 
cosine  factor  approaches  zero,  since: 

lim  sin  x  _  ^ 

X  — >  co  X 

and  the  intensity  becomes  equal  to  A§/4  everywhere  except  at  the  position  of 
the  maximum  of  the  envelope.  At  this  point  the  argument  of  the  sine  is  zero, 
regardless  of  the  value  of  (vg  -  v^) ,  and  the  former  limit  given  in  Eq.  45 
holds.  Thus,  at  the  position  of  the  maximum  of  the  envelope,  the  maximum 
intensity  possible  is  A§/2,  subject  to  modulation  by  the  cosine  factor,  and  ' 
only  one  fringe  remains. 

Most  important  however  is  the  form  of  the  indices  of  refract¬ 
ion  which  appear  in  Eq.  44.  Assuming  for  the  moment  that  no  electrons  are 
present,  (the  effect  of  electrons  was  considered  earlier)  it  is  noted  that 
the  shift  of  individual  fringes  with  frequency  v  is  controlled  by  AQ  and 
(l  +  A' ),  while  the  shift  of  the  envelope  with  frequency  (v2  - ■ v, )  is  con¬ 
trolled  by  (Aq  -  A-^)  and  (l  +  A^  +  A£).  Thus,  as  mentioned  earlier,  the 
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envelope  will  be  shifted  by  a  different  amount  than  will  be  individual 
fringes  for  a  given  change  in  (p2  *  P^).  This  effect  is  known  as  "rollover" 
(Ref.  5)>  and  the  reason  for  its  occurrence  is  shown  in  Appendix  A. 


3.  EXPERIMENTAL  RESULTS  FOR  OXYGEN 
l)  Procedure  and  Method  of  Analysis 

The  basic  method  which  was  used  to  determine  the  Gladstone- 
Dale  constants  has  been  outlined  in  Appendix  B  of  Ref.  1.  This  involved 
passing  a  strong  shock  wave  through  the  gas  under  study  and  measuring  the 
index  of  refraction  of  the  hot  dissociating  gas  behind  the  shock  wave.  This 
could  be  done  in  the  UTIAS  4"  x  7"  Shock  Tube,  using  the  9"  Mach-Zehnder 
interferometer  to  measure  the  index  of  refraction.  A  convenient  reference 
index  of  refraction  was  that  of  the  cold  gas  ahead  of  the  shock  wave,  for 
which  all  the  parameters  of  state  were  known,  while  the  change  in  the  index 
of  refraction  between  the  shocked  and  unshocked  gas  could  be  determined  from 
an  instantaneous  interferogram  of  the  shock  wave  as  it  passed  through  the 
field  of  the  interferometer. 


As  noted  above,  all  the  required  parameters  of  the  gas 
ahead  of  the  shock  wave  were  known.  Also,  the  Gladstone-Dale  constants  at 
various  wavelengths  for  various  gases  have  been  determined  from  room- 
temperature  conditions,  and  were  available  in  Ref.  6.  For  the  shock- 
heated  gas,  the  required  state  parameters  could  not  be  determined  directly. 
Consequently,  values  of  the  density  ratio  across  the  shock  wave,  had  to  be 
determined  from  the  initial  conditions,  using  normal  shock  calculations. 
Using  this  data,  it  was  then  possible  to  determine  the  Gladstone-Dale  con¬ 
stants  for  the  hot  gas  from  the  refractive  index  measurements. 


For  the  case  of  oxygen,  the  contribution  of  the  electrons 

produced  by  the 
neglected  (  u)|  / 
for  the  shocked  gas  becomes: 


shock  wave,  to  the  total  refractivity  of  the  gas,  can  be 
2u>2  r:  10-16) ,  s0  ^at  the  Gladstone-Dale  equation  (Eq.  34) 
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i  =  n,p. 
i  1  i 


or 

n2  -  1  -  p2  [%  (1  -  a.)  +  Ka  OgJ  (45) 


where  p2  is  the  total  density  of  the  shocked  gas,  a2  is  the  degree  of 
dissociation,  i.e.,  the  fraction  of  the  oxygen  molecules  that  are  split 
up  into  their  component  atoms,  %  is  the  Gladstone-Dale  constant  for  the 
oxygen  molecule,  and  Ka  is  the  Gladstone-Dale  constant  for  the  oxygen  atom. 

Similarly,  for  the  cold  gas  ahead  of  the  shock  wave,  Eq.  34 

becomes: 

n1-l=p1  ^  (46) 


where  Ko2  is  the  Gladstone-Dale  constant  of  the  oxygen  molecule  at  room 
temperature,  as  given  in  Ref.  6.  On  the  basis  of  the  arguments  present  in 
Sec.  2.1.5,  Kq2  would  be  expected  to  be  equal  to  %,  in  spite  of  the  vastly 


different  temperatures  involved  in  the  two  determinations.  However,  one  of 
the  objectives  of  this  experimental  work  is  to  check  this  point,  and  there¬ 
fore,  for  the  sake  of  generality,  no  assumption  is  made  regarding  the  value 
of  %. 


The  relation  between  the  non-dimensional  fringe  shift  and 
the  change  in  index  of  refraction  is  given  by  (Eq.  4o): 

TT  ^nP  "  V  =  \  L  (yP  "  Vq)  =  d  (yP  "  yQ^  =  SQP  ^ 

If  nq  represents  the  initial  state  1  ahead  of  the  shock  wave  on  the  interfero- 
gram,  and  np  represents  the  final  state  2  behind  the  shock  wave,  one  can  sub¬ 
stitute  from  Eqs.  45  and  46  for  np  and  n^  respectively  and  obtain: 

r  p2  k  (1  -  <%)  +  h  “si  -  T  P1  =  S12 

o  •-  Jo 

or 

%[^p2(i.a2l  ♦  %  fr  ps  “s]  =  S12  +  pi  ko2 
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Equation  47  contains  two  unknowns,  so  that  one  determination 
of  SjL2  would  not  provide  enough  information.  A  minimum  of  two  equations 
of  this  form  are  required  to  determine  both  %  and  KA,  i.e.,  two  independent 
runs  in  the  shock  tube  are  required. 

This  procedure  of  using  two  simultaneous  equations  to  deter¬ 
mine  %  and  KA  was  that  suggested  in  Appendix  B  of  Ref.  1.  However,  it  had 
to  be  discarded  since  it  was  found  that  the  values  obtained  for  Kj^  and  KA 
were  extremely  sensitive  to  small  errors  in  the  measurement  of  the  fringe 
shifts.  An  alternative  procedure  was  that  used  in  Ref.  7,  where  was 
arbitrarily  set  equal  to  Kq2»  thus  eliminating  one  unknown.  This  allowed 
Ka  to  be  determined  from  only  one  measurement  of  fringe  shift,  and  the 
error  in  KA  resulting  from  an  error  in  S^2  was  much  reduced.  However, 
while  on  the  basis  of  the  classical  arguments  presented  earlier,  there  is 
some  Justification  for  setting  equal  to  Kq_,  until  this  equality  is 
actually  proved  to  be  so,  the  value  obtained  for  K&  is  still  open  to  question, 
since  it  is  influenced  by  the  choice  of  the  value  of  K$j. 

Obviously,  it  would  be  desirable  to  remove  this  assumption, 
and  so  a  new  procedure*  was  developed.  It  was  better  in  that  it  removed 
the  assumption  of  a  value  for  %,  and  allowed  a  large  number  of  determinations 
of  S]_2  to  be  used  to  obtain  one  pair  of  values  for  %  and  K^.  This  method 
is  outlined  below. 


*  The  author  is  indebted  to  Dr.  A.  Levy  for  the  original  suggestion  leading 
to  the  development  of  this  method. 


Equation  47  can  be  written  in  the  form: 


r;  [%  - a2>  *  ka “a]  - sig  *“iTo\  m 

r  p2  lK'i  -  sia  *  ’i  r  “oa 

O  0 

K’  =  %  (1  -  a2)  +  KA  a2  (48a) 

It  is  seen  that  K'  is  the  specific  refractivity  of  the  mix¬ 
ture  of  O2  and  0,  and  the  value  of  K*  changes  depending  on  the  value  of  0^. 

In  particular,  at  C*2  =  K'  becomes  equal  to  K^,  and  at  02  =  K'  becomes 

equal  to  KA.  Also,  if  Eq.  48a  is  written  in  the  form: 

K’  =  %  +  a2  (Ka  -  %)  (48b) 

it  is  seen  that  the  value  of  K'  varies  linearly  with  ou,  since  KA  and  are 
constants  for  the  gas  under  study.  Thus,  if  K’  is  plotted  versus  Qf2>  anc*  a 
straight  line  is  drawn  through  the  points,  it  is  seen  that  the  intercept  at 
(*2  =  0  gives  a  value  of  and  the  intercept  at  a2  =  1  gives  a  value  of  KA« 
The  value  of  this  procedure  is  that  there  is  no  limit  to  the  number  of 
points  K*  versus  a.2  that  can  be  used,  i.e.,  the  number  of  runs  in  the  shock 
tube,  since  each  run  produces  one  point  of  K’  versus  a„,  and  every  additional 
point  or  run  reduces  the  uncertainty  in  the  slope  and  intercepts  of  the  line 
drawn  through  the  points. 

•  The  criteria  for  choosing  running  conditions  in  the  shock 
tube  can  be  obtained  with  the  aid  of  Eq.  48.  Solving  for  K'  in  Eq.  48,  one 
obtains : 


and  if  one  assumes  that  pj  ,  L,  >Q,  and  Kq2  are  known  exactly,  one  can  write: 


o  o 


where  AK’,  AS12,  anti  Ap2  are  errors  in  these  quantities,  and  where 
AK'  is  the  sum  of  the  errors  in  the  terms  involving  AS12  and  AP2  •  Thus 
it  is  seen  that  the  error  in  K*  may  be  reduced  by  raising  the  density  of  the 
shocked  gas,  Pg.  Also,  since  it  is  the  intercepts  at  0-2  =  0  and  a2  =  1 
that  are  wanted,  it  is  desirable  to  obtain  points  with  02  as  low  as  possible 
and  also  with  OC2  as  high  as  possible.  Obviously,  the  limits  on  the  choice 
of  running  conditions  are  set  by  the  structural  strength  of  the  shock  tube. 
In  addition,  it  was  desired  to  keep  the  temperature  of  the  shocked  gas  with¬ 
in  the  range  of  approximately  3300°K  to  4400°K,  which  was  roughly  the  tem¬ 
perature  range  involved  in  the  oxygen  expansion  studies  (Ref.  l)  which  the 
present  work  supports.  The  change  with  temperature  of  the  Gladstone-Dale 


or: 

where : 
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constants  would,  then  be  shown  by  the  difference  between  fQ-gj  the  rcom- 
temperature  value  for  the  oxygen  molecule,  and  %,  the  value  obtained  in  the 
present  work  at  the  elevated  temperature.  Also,  any  variation  of  the 
Gladstone-Dale  constants  with  temperature  within  the  range  chosen  for  this 
work  would  show  up  as  scatter  in  the  points  K*  vs  a?,  and  the  scatter  would 
show  some. sort  of  trend  with  temperature. 

The  best  way  to  fit  a  straight  line  to  the  values  of  K'  vs 
£*2  is  to  use  a  least  squares  analysis;  however,  the  procedure  is  not  as 
straight  forward  as  usual  in  this  case.  Normally,  in  minimizing  the  sum  of 
the  squares  of  the  deviations,  one  assumes  that  the  variances  of  all  the 
points  are  equal.  This  would  be  unrealistic  in  the  present  case,  since  the 
variances  are  in  fact  not  equal,  as  indicated  by  the  estimated  error  (See 
Table  3)*  Consequently,  a  quantity  p  is  calculated  for  each  point  such  that 
p  is  proportional  to  the  inverse  of  the  estimated  percent  error  of  the  poinjj. 
Thus,  instead  of  one  point  at  each  location,  one  assumes  p2  points,  where  p^ 
is  not  necessarily  an  integer.  The  equations  giving  the  slope  (KA  -  %)  and 
the  intercept  %  at  ~  0  then  become  (see  Eq,  ^8b) 

(ka  -  Kj  £  p2  a2  -  £p2  a  K*  +■  p2  a  =  o 

(50) 

(Ka  -  %)  £  p2  a  -  £p2  K'  +  %  £p2  =  0 
These  equations  could  al$o  be  solved  for  KA  and  explicitly. 


obtained  from: 


p 

where:  pa 


where  6  is  the  percent  error  of  the  point  at  abscissa  a  and  5'  is  the  per¬ 
cent  error  of  the  point  (possibly  hypothetical)  for  which  the  quantity  p  is 
equal  to  one.  The  expression  in  Eq.  51  was  originally  derived  (Ref.  15)  for 
standard  deviation,  but  it  is  assumed  here  that  the  percent  error  is  propor¬ 
tional  to  the  standard  deviation  and  that  they  remain  porportional  throughout 
the  calculations , 

(Equations  50  and  51  are  actually  adaptations  of  expressions 
derived  for  the  case  of  equal  variance  in  all  the  points.  However,  once  the 
method  for  weighting  the  points  has  been  decided,  the  forms  shown  in  Eqs.  50 
and  51  follow  immediately.  The  particular  method  of  weighting  used  here  was 
decided  upon  so  that  the  results  obtained  for  6  using  Eq.  51  would  be  inde¬ 
pendent  of  the  actual  numerical  magnitude  of  p,  the  only  requirement  being 
that  p  be  proportional  to  the  inverse  of  the  estimated  percent  error,  as  de¬ 
tailed  earlier.  The  condition  that  £p2  must  be  equal  to  the  number  of 
points  is  required  for  the  test  of  significance  which  is  applied  in  Sec.  1.2 
to  determine  whether  the  difference  between  the  results  for  %  and  KA  is 
significant  in  view  of  the  experimental  error  (See  Ref.  15).  The  result  of 
this  test  is  influenced  by  the  magnitude  of  £p2  and,  since  in  the  original 
equal- variance  version  of  this  test  the  quantity  corresponding  to  Ep2  was 
the  number  of  experimental  points,  it  was  decided  to  set  £p2  equal  to  the 
number  of  points). 


A  measure  of  the  percent  error  at  the  intercepts  can  be 
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3.2  Experimental  Results  and  Error  Analysis 

A  summary  of  the  runs  in  the  UTIAS  4"  x  7"  Shock  Tube  used 
for  the  determination  of  the  Gladstone-Dale  constants  for  dissociating  oxygen 
is  shown  in  Table  1.  Twelve  runs  were  taken:  shock  Mach  numbers  varied  be¬ 
tween  8.86  and  16.64,  while  initial  pressure  varied  from  3*99  mm  Hg  to 
54.0  mm  Hg.  It  is  also  noted  that  the  initial  temperature  did  not  vary  more 
than  one  degree  Kelvin  from  300°K. 

As  mentioned  before,  the  state  parameters  of  the  gas  behind 
the  shock  vave  had  to  be  determined  from  normal  shock  calculations,  and  it 
was  plain  that  any  results  obtained  for  the  Gladstone-Dale  constants  would 
be  only  as  accurate  as  the  calculations  used.  Fortunately,  suitable  normal 
shock  data  was  obtainable  from  Ref.  13.  This  data  was  calculated  for  T-,  of 
300°K,  which  was  in  line  with  the  actual  T]_  in  the  shock  tube.  From  this 
data,  values  of  the  degree  of  dissociation  of  the  shock-heated  gas,  the  den¬ 
sity  ratio  across  the  shock  wave,  and  the  equilibrium  temperature  of  the  hot 
gas  were  obtained  for  each  run. 

Of  major  importance  too  was  the  accuracy  with  which  measure¬ 
ments  of  fringe  shift  on  the  interfercgrams  could  be  mede.  Prints  made 
irora  typical  inter ferograms  are  shown  in  Figs.  1,  2,  and  3*  Since  the 
fringes  were  discontinuous  across  the  shock  wave,  use  had  to  be  made  of  a 
white-light  interferogram,  taken  at  the  same  time  as  the  monochromatic 
interferogram,  in  order  to  trace  the  zero-order  fringe  across  the  shock 
front.  This  was  possible  because,  as  shown  in  Figs.  1  and  2,  only  a  few 
fringes  on  either  side  of  the  zero-order  fringe  are  visible  on  a  white  light 
interferogram.  However,  the  phenomenon  of  "rollover  ",  which  becomes 
important  when  there  are  large  density  changes,  meant  that  great  care  had  to 
be  exercised  when  using  the  white-light  fringes  to  identify  the  zero-order 
fringe  (See  Sec.  2.2.4  and  Appendix  A). 

In  addition,  care  had  to  be  exercised  when  measuring  the 
fringe-shift  since  the  shock  front  was  so  badly  blurred  on  the  earlier  inter- 
ferograms  such  as  in  Fig.  1.  These  interferograras  were  taken  with  a 
magnesium  spark  source  in  the  interferometer,  which  limited  the  time  resolu¬ 
tion  on  the  interferograms  to  about  two  microseconds.  Consequently,  account 
had  to  be  taken  of  the  curvature  of  the  fringes  between  the  two  points  of 
measurement  in  front  of  and  behind  the  shock  wave,  since  these  points  were 
far  apart  due  to  the  apparent  thickness  of  the  shock  front.  Later  inter¬ 
ferograms  were  taken  with  the  exploding  wire  light  source  and  a  Kerr  Cell 
(See  Appendix  C)  which  improved  the  time  resolution  by  a  factor  of  10.  The 
difference  in  sharpness  of  the  shock  front  is  readily  apparent  in  the  inter¬ 
ferograms  shown  in  Fig.  3* 

Another  difficulty  which  had  to  be  contended  with  is  appar¬ 
ent  in  Fig.  lb.  This  was  the  so-called  "phasing"  of  the  white-light  fringes, 
in  which  the  fringe  envelope  of  the  white-light  fringes  moved  with  respect 
to  individual  fringes  as  one  followed  the  fringes  from  one  side  of  the 
interferogram  to  the  other.  This  phasing  was  very  similar  to  rollover  ex¬ 
cept  that  it  was  a  function  of  the  optics  of  the  interferometer,  not  of  a 
density  change,  since  phasing  was  present  on  the  no-flow  as  well  as  the  flow 
interferograms.  As  a  result,  determination  of  the  position  of  the  zero- 
order  fringe  from  an  estimate  of  the  rollover  on  the  white-light  interfero¬ 
grams  was  a  very  difficult  matter.  Eventually,  the  problem  was  overcome  by 


taking  a  series  of  densitometer  traces  across  the  white-light  fringe  pattern 
at  various  positions  on  both  the  flow  and  no- flow  interferograms .  Since 
the  phasing  was  present  on  both  the  flow  and  no-flow  interferograms,  while 
the  rollover  was  present  only  on  the  flow  interf erograms,  it  was  possible 
to  separate  the  "rollover"  from  the"phasing"  using  the  densitometer  traces, 
and  thus  identify  the  zero-order  fringe. 

No  explanation  for  this  phasing  could  be  advanced.  Eventually 
it “disappeared  and  was  replaced  by  another  effect  which  is  apparent  in  Fig. 

2.  Here,  the  spacing  of  the  fringes  varied  across  the  white-light  pattern, 
and  the  fringe  envelope  was  not  symmetrical  about  its  maximum.  Although 
the  interferograms  shown  in  Fig.  2  were  taken  with  the  exploding  wire  source 
and  Kerr  Cell,  the  disappearance  of  the  phasing  and  the  appearance  of  the 
asymmetric  white-light  fringe  pattern  did  not  correspond  to  the  introduction 
of  the  exploding  wire  light  source;  this  change  was  apparent  in  other  work 
done  with  the  magnesium  spark  source  before  the  exploding  wire  light  source 
was  introduced.  Again,  no  explanation  could  be  advanced  for  this  effect. 

However,  the  disappearance  of  the  phasing  made  it  possible 
to  determine  the  amount  of  rollover  by  inspection  since  it  was  no  longer 
complicated  by  other  effects,  and  thus  identification  of  the  zero-order 
fringe  became  relatively  simple.  This  was  fortunate,  because  these  particu¬ 
lar  interferograms  were  taken  with  the  exploding  wire  source,  and  grain  pro¬ 
blems  resulting  from  the  use  of  extremely  high-speed  film  (See  Appendix  C), 
along  with  the  fact  that  so  many  fringes  were  visible  in  the  white-light 
pattern,  made  densitometer  traces  taken  on  these  interferograms  rather  hard 
to  interpret. 


After  the  position  of  the  zero-order  fringe  had  been  deter¬ 
mined  on  both  sides  of  the  shock  wave  for  the  white-light  interferograms, 
the  dimensional  fringe  shift  was  determined  using  a  travelling  microscope. 
This  measurement  was  then  transferred  to  the  monochromatic  interferogram 
and  accurate  measurements  of  the  dimensional  fringe  shift  and  of  the  fringe 
spacing  were  made,  from  which  the  non-dimensional  fringe  shift',  S-^,  was 
determined.  It  was  found  that  the  dimensional  fringe  shifts  on  tne  white- 
light  and  monochromatic  interferograms  were  never  quite  the  same.  This 
caused  some  difficulty  at  first  until  it  was  realized  that,  unless  the  wave¬ 
length  of  the  monochromatic  interferogram  and  the  wavelength  corresponding 
to  the  spacing  of  individual  fringes  on  the  white-light  interferogram  were 
the  same,  dispersive  effects  would  result  in  the  dimensional  fringe  shifts 
being  slightly  different.  In  addition,  the  monochromatic  and  white-light 
interferograms  were  taken  with  different  camera  backs,  so  that  the  magnifi¬ 
cations  were  likely  to  be  slightly  different,  due  to  small  differences  in 
focussing.  It  ras  also  found  to  be  of  some  importance  to  identify  the  zero- 
order  fringes.  This -was  necessary  because  it  was  found  that  the  monochro¬ 
matic  fringe  spacing  tended  to  change  slightly  across  the  field,  possibly 
due  to  imperfections  in  the  interferometer  optics,  together  with  small 
vibrations  transferred  to  the  interferometer  from  the  building.  This  change 
in  fringe  spacing  was  not  important  in  itself,  but  it  meant  that,  unless  one 
transferred  to  dimensional  fringe  shift  on  the  white-light  interferograms 
to  the  corresponding  fringes  on  the  monochromatic  interferograms,  the 
dimensional  white-light  fringe  shift  would  not  correspond  to  actual  fringes 
on  the  monochromatic  interferograms,  leading  to  possible  errors. 


27 


When  all  factors  were  taken  into  account,  it  was  judged  that 
the  error  in  the  fringe  shift  measurements  was  about  0.1  fringe.  Measure¬ 
ments  on  the  earlier  interferograms  taken  with  the  magnesium  spark  source 
could  be  repeated  to  better  than  0.1  fringe,  but  it  was  felt  that  allowance 
had  to  be  made  for  systematic  errors  arising  from  the  blurring  of  the  shock 
front.  For  the  interferograms  taken  with  the  exploding  wire  source  and  the 
Kerr  Cell,  it  was  judged  that  systematic  errors  were  much  reduced  due  to  the 
enhanced  crispness  of  the  fringes  and  the  better  definition  of  the  shock 
front.  However,  due  to  variations  in  film  and  in  processing,  some  of  the 
interferograms  were  not  of  the  best  quality,  and  overall  accuracy  of  better 
than  0.1  fringe  could  not  be  claimed. 

Having  deterimed  the  fringe  shifts,  it  was  then  possible  to  cal¬ 
culate  values  of  K'  and  the  details  of  these  calculations  are  shown  in  Table 
2.  The  values  of  the  parameters  which  entered  into  the  calculations  are  also 
shown.  In  particular,  the  initial  density,  p^,  was  calculated  from  the 
pressure  and  the  temperature  using  the  perfect  gas  law,  a  good  approximation 
at  the  low  temperatures  and  pressures  involved.  The  values  of  K'  were  then 
plotted  as  shown  in  Fig.  4. 


Before  a  least  squares  analysis  could  be  carried  out  to  find  the 
best  straight  line  through  the  points,  an  estimate  had  to  be  made  of  the  error 
to  be  expected  in  K' .  How  these  errors  enter  is  shown  in  Eq.  49.  As  mention¬ 
ed  above,  the  error  in  S12  was  estimated  as  0.1  fringe.  A  further  error 
entered  in  due  to  uncertainties  in  the  measurement  of  shock  speed,  and  the 
effects  of  this  are  shown  in  Table  3»  where  estimates  h\ve  been  made  of  the 
resultant  errors  in  a2  and  Pp/Pl*  Fortunately,  a2  does  not  enter  into  the 
calculation  of  K',  but  it  does  determine  the  position  of  the  K'  vs  a2  plot 
used  to  calculate  the  intercepts.  However  it  is  seen  from  Fig.  4  that  the 
error  in  Og,  ranging  from  .004  to  .015  is  small  compared  to  the  error  in  K' 
itself. 


It  is  seen,  however,  that  the  errors  in  K'  resulting  from  errors 
in  P2/P1  are  comparable  to  the  errors  resulting  from  uncertainties  in  S^2 
and  points  up  the  necessity  of  accurately  determining  shock  speed.  The  de¬ 
tails  of  these  calculations  are  shown  in  Table  3» 


I 
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It  was  determined  that  errors  in  K’  and  ctp  arising  from  uncer¬ 
tainties  in  the  initial  pressure  could  be  neglected.  The  uncertainty  in  p^ 
was  -.05  mm  Hg. 

Having  then  determined  the  total  error  in  K‘ ,  which  ranged  from 
.80#  to  2.78$,  a  least  squares  analysis  was  carried  out  using  Eqs.  50,  and 
the  result  of  the  calculation  is  shown  in  Fig.  4.  It  is  noted  that  the 
deviations  of  the  points  from  the  line  agree  very  well  with  the  errors  calcu¬ 
lated  on  the  basis  of  estimated  error  in  fringe  shift  and  shock  speed  mea¬ 
surements.  The  actual  deviations  of  the  points  K’  from  the  line  are  listed 
in  Table  3  for  comparison  with  the  estimated  errors. 

Finally,  the  errors  in  the  intercepts  were  determined  using  Eq. 
51,  leading  to  the  final  results: 

I'm  *  i-»  x  if4  -if- 

Ka  =  2.04  X  10-U  — -  ±  2$ 
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The  value  here  obtained  1'or  Kj^  can  be  compared  to  1.902  x  10 
m3/kg,  the  value  of  the  Gladstone-Dale  constant  for  oxygen  at  room  temperature 
as  given  in  Ref.  6.  It  is  noted  that  the  value  given  in  Ref.  6  falls  only  just 
outside  the  experimental  error  in  Ihe  present  result.  This  agreement  was  con¬ 
sidered  to  be  quite  good,  considering  the  radically  different  methods  used  to 
obtain  the  two  values.  This  agreement  would  also  tend  to  confirm  the  conclusion 
reached  in  Sec.  2.1.5  that  the  Gladstone-Dale  constants  for  oxygen  do  not  change 
with  temperature.  Greater  certainty  in  this  respect  will  have  to  await  new 
developments  leading  to  reduction  in  the  errors  in  and  P2/Pi* 


3  /  + 

The  value  obtained  for  KA  can  be  compared  to  1.8  x  10  m  /kg  - 
12$  as  obtained  in  Ref.  7.  However,  the  work  of  Ref.  7  was  carried  out  using  a 
smaller  shock  tube  and  much  lower  densities  than  were  used  in  the  present  work, 
so  that  the  experimental  error  was  much  higher  than  in  the  present  work.  Also, 
as  mentioned  earlier,  the  result  of  Ref.  7  for  Ka  assumed  that  Km  was  equal  to 
*02’  thereby  casting  doubt  on  the  reliability  of  the  value  obtained  for  K^.  In 
view  of  the  present  work,  however,  it  would  appear  that  the  assumption  of  % 
equal  to  Ko2  was  valid,  so  that  the  present  result  for  and  that  of  Ref.  7  can 
be  considered  on  roughly  equal  footing.  It  is  then  noted  that,  when  the  experi¬ 
mental  errors  are  taken  into  account,  the  two  results  sire  in  agreement.  The 
result  of  the  present  work  however  is  considered  more  reliable  in  view  of  the 
much  lower  experimental  error. 


It  was  noted  that  the  present  work  indicated  that  the  Gladstone- 
Dale  constant  for  the  atom  was  greater  than  that  for  the  molecule,  rather  than 
smaller  as  indicated  by  the  result  of  Ref.  7-  As  a  further  check  on  this,  a 
statistical  test  was  applied  to  the  graph  shown  in  Fig.  4  to  determine  whether 
the  slope  of  K'  vs  a2  was  significant  in  view  of  the  scatter  of  the  points  and 
the  experimental  error  in  each  point  (see  Ref.  15).  The  hypothesis  tested  was 
that  the  slope  K'  versus  a2  was  zero,  and  it  was  found  that,  at  5$  significance 
level,  the  slope  of  the  line  K’  versus  Og  was  significant*,  indicating  that  KA 
was  indeed  greater  than  Km* 


As  a  further  check  on  the  validity  of  the  results,  two  tests 
were  applied  to  ascertain  that  there  were  indeed  no  temperature  effects  and  no 
density  effects.  These  effects  (if  any)  were  ignored  in  the  analysis  used  to 
obtain  %  and  KA,  but  if  they  were  present,  their  effect  would  show  up  as  scatter 
in  the  points,  and  the  scatter  would  be  correlated  in  some  way  with  the  tempera¬ 
ture  and/or  density.  Therefore,  the  deviation  of  each  point  from  the  mean  line 
in  Fig.  4  has  been  plotted  against  temperature  T2  in  Fig.  5,  and  against  density 
in  Fig.  6.  The  scatter  appears  to  be  completely  random  in  both  cases,  indicating 
that  temperature  and  density  effects  are  not  present. 


*  This  means  that,  if  the  slope  were  indeed  zero,  and  the  twelve  runs  were  re¬ 
peated  an  infinite  number  of  times,  the  slope  of  the  line  K'  versus  (*2  would 
assume  many  different  values,  averaging  to  zero.  The  statistical  test  arrives 
at  two  extreme  values  for  the  slope  such  that  95$  of  all  the  values  obtained  for 
the  slope  would  fall  within  the  range  defined  by  the  two  extremes.  In  the  pre¬ 
sent  case,  the  actual  value  of  the  slope  falls  outside  this  range,  and  the  con¬ 
clusion  is  reached  that  the  slope  is  significant,  i.e.,  that  the  slope  is  differ¬ 
ent  from  zero  . 
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4.  CONCLUSIONS 

An  experimental  determination  of  the  Giadstone-Dale  constants 
for  molecular  and  atomic  oxygen  was  carried  out  using,  the  UTIAS  4"  x  7"  Hyper- 
sonc  Shock  Tube.  The  values  obtained  were  1.93  x  10“4  m3/kg  t  1 $  for  molecular 
oxygen,  and  2.04  x  ICT4  m3/kg  -  2 $  for  atomic  oxygen. 

The  result  obtained  for  the  Gladstone-Dale  constant  of  the  oxygen 
molecule  was  in  good  agreement  with  the  accepted  room  temperature  value  (Ref.  6) 
of  1.902  x  10”4  m^/kg.  This  appeared  to  confirm  the  conclusion,  based  on 
classical  arguments,  that  the  Gladstone-Dale  constants  for  oxygen  do  not  change 
with  temperature . 

The  value  previously  obtained  by  Alpher  and  White  (Ref.  7)  for 
the  Gladstone-Dale  constant  of  the  oxygen  atom  was  1.6  x  ICT4  m3/kg  -  12$. 

However,  the  wuik  of  Ref.  7  involved  the  assumption  of  a  value  for  the  specific 
refractivity  of  the  oxygen  molecule.  This  condition  was  removed  in  the  present 
work,  so  that  the  value  here  obtained  for  the  specific  refractivity  was  free  of 
any  assumptions.  In  addition,  the  experimental  error  in  the  present  work  was 
much  lower  than  in  Ref.  7* 

A  number  of  tests  were  carried  out  to  test  the  validity  of  the 
results  obtained.  A  statistical  test  was  carried  out  to  determine  whether  the 
difference  in  the  results  obtained  for  the  specific  refractivities  of  the  molecule 
and  «.tom  was  significant  in  view  of  the  experimental  error  in  the  points;  the 
result  showed  that  the  refractivities  were  indeed  different.  Tests  applied  to 
determine  whether  there  was  any  temperature  or  density  effects  in  the  results 
showed  that  these  effects  were  absent. 

In  addition  to  the  experimental  work,  a  theoretical  analysis  was 
carried  out  in  which  expression  for  the  index  of  refraction  of  dissociating  and 
ionizing  gases  were  derived  from  Maxwell's  equations  and  the  equations  of  motion 
for  a  plasma.  The  role  of  the  Gladstone-Dale  constant  was  shown  along  with  its 
relation  to  the  polarizability  of  molecules  and  atoms.  A  further  analysis  was 
then  done  to  show  the  application  of  the  indices  of  refraction  to  interferometry 
and  to  show  analytically  how  monochromatic  and  white-light  interference  patterns 
are  formed. 
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TABLE  1 


Summary  of  Runs  In  the  UIIAS  4"  x  7"  Shock  Tube  Used  to  Determine  the  Gladstone- 

Pale  Constants 


Run 

Ms 

pl 

T1 

a2 

p2/pi 

T2 

114 

13.48 

8.05 

mm  Hg 

299. 6°K 

.343 

13.41 

3935°K 

115 

12.46 

17.4 

t» 

300.4  " 

.256 

12.21 

3925  " 

116 

11.61 

24.2 

II 

299.8  " 

.195 

11.34 

3835  " 

117 

11.99 

15.0 

II 

299.4  " 

.226 

11.90 

3810  " 

118 

16.64 

3.99 

II 

299.3  " 

.630 

15.58 

4295  " 

119 

15.12 

■  9.94 

II 

300.0  " 

.473 

14.27 

4265  " 

120 

9.57 

38.0 

It 

299.1  " 

.0793 

9.27 

3460  " 

121 

10.19 

12.6 

It 

300.1  " 

.122 

10.27 

3440  " 

197 

8.86 

54.0 

It 

300.9  " 

.0496 

8.44 

3310  " 

198 

13.20 

21.5 

II 

299.3  " 

.304 

12.63 

4no  " 

200 

15.75 

7.50 

II 

300.6  " 

.534 

14.77 

4300  " 

201 

15.64 

7.50 

II 

299.9  " 

.524 

14.72 

4280  ” 

TABLE  2 


Calculation  of  K' 


Run 

P1 

P2 

3 12 

K* 

114 

1.378  X 

10-2 

kg/m^ 

.1848  kg/rn^ 

6.50 

1.945  X 

-4 

10 

m^/kg 

115 

2.971 

ft 

.3628 

ft 

13.00 

1.993 

tf 

ft 

116 

4.l4o 

ft 

.4695 

ft 

16.50 

1.970 

If 

ft 

117 

2.570 

If 

.3058 

It 

10.88 

1.984 

tf 

ft 

118 

.6838 

tf 

.1065 

ft 

3.78 

1.942 

ft 

ft 

119 

1.699 

ft 

.2424 

ft 

8.70 

1.974 

It 

tf 

120 

6.516 

ft 

»• 

.6o4o 

ft 

20.20 

1.920 

ft 

ft 

121 

2.153 

ft 

.2211 

ft 

7.53 

1.931 

ft 

ft 

197 

9.205 

!» 

.7769 

ft 

26.00 

1.942 

ft 

tf 

198 

3.648 

ft 

.4653 

ft 

16.40 

1.958 

tf 

ft 

200 

1.280 

ft 

.1891 

ft 

6.83 

1.981 

tf 

ft 

201 

1.283 

If 

.1889 

It 

6.94 

2.013 

tf 

ft 

TABLE  3 


Error  in  Values  of  K1 


Run 

Error  in  M 

s 

absolute 

percent 

Error  in  a2 

Error  in  p2/ 

114 

±  .10 

±  .73 

±  .008 

±  .069 

115 

.08 

.68 

.006 

.065 

116 

.07 

.63 

.004 

.061 

117 

.08 

.65 

.006 

.067 

118 

.15 

•  90 

.015 

.047 

119 

.12 

.82 

•  .011 

.060 

120 

.05 

.52 

.002 

00 

-3- 

0 

• 

121 

.06 

.55 

.003 

.062 

197 

.04 

.48 

.002 

.039 

198 

.09 

.71 

.007 

.059 

200 

.13 

.85 

.012 

.059 

201 

.13 

.85 

.012 

.059 

ffc 


TABLE  3  (con'd) 


Run 

A  K' 

AK'/K' 

observed  AK' 
(from  Fig. 4) 

114 

±3.78  X 

10  ^  m3/kg 

1.94 1o 

2.50  x  10"6  m3 

115 

2.47 

tl  It 

1.24  56 

3.20 

116 

2.15 

It  l» 

1.0956 

1.52 

117 

2.79 

II  II 

l.4l" 

2.60  " 

118 

5.40 

II  II 

2.7 8$ 

5.74 

119 

2.95 

It  II 

1.49 56 

0.93 

120 

GO 

l»  II 

.95856 

2.29  " 

121 

3.49 

II  It 

1.8l% 

1.63 

197 

1.56 

II  II 

.80356 

.21 

198 

2.02 

II  II 

1.0356 

.80 

200 

3.50 

II  II 

1.77% 

.85 

201 

3.52 

It  II 

1.7556 

2.45 

i 


FIG.  I  -  TYPICAL  INTERFEROGRAMS  TAKEN  WITH  MAGNESIUM  SPARK  LIGHT  SOURCE  -  SHOWN  APPROXIMATELY  3/4  ACTUAL  SIZE 

(Run  H9:  M  *  15.12,  p  *  9.94  mm.Hg.  Oxygen) 


-^4-  i 


(Run  200;  M  *  15.  75,  p  «  7.50  mtn.  Hg.  Oxygen) 


FIG.  3a  -  ENLARGEMENT  OF  MONOCHROMATIC  FLOW  INTERFEROGRAM 
SHOWN  IN  FIG.  1  (Approximately  twice  full  size) 


FIG.  3b  -  ENLARGEMENT  OF  MONOCHROMATIC  FLOW  INTERFEROGRAM 
SHOWN  IN  FIG.  2  (Approximately  four  times  full  s  ize) 


FIG.  4  PLOT  OF  K*  vs M  USED  TO  CALCULATE  VALUES  OF  GLADSTONE-DALE  CONSTANTS 


‘PARABOLIC  MIRROR  I 


\ 


MIRROR  2 


FIG.  7  -  GEOMETRY  CF  9-INCH  MACH-ZEHNDEK  INTERFEROMETER  (TAKEN  FROM  REF.  3) 


FIG.  8  -  SHIFT  OF  PARALLEL  FRINGES 


APPENDIX  A 


Calculation  of  White-Light  Interference  Patterns 

In  this  appendix,  details  of  the  integration  of  Eq.  43  leading 
to  Eq.  44  will  be  given.  In  addition,  the  integration  will  be  carried  out  for 
the  case  of  an  intensity  distribution  which  is  weighted  toward  one  end  of  the 
frequency  spectrum,  rather  than  the.  symmetrical  distribution  assumed  in  Eq. 

44.  Details  will  be  given  of  numerical  calculations  carried  out  to  illustrate 
the  effects  noted  in  the  calculations. 


The  intensity  distribution  in  a  white-light  interference  pattern 
is  given  basically  by  Eq.  M; 


2  — 

I  dv 


Aq 


I 


2  2; 
cos 


7TLvK 

c 


<V  h1  ■ 


^  27reyvna 
2  c 


dv  (43) 


Eor  simplicity,  we  substitute: 

Y  -  ^4  .  Y  - 

“1  c  ’  A2  ~  c 

As  noted  in  Sec.  2.2.4.,.  the  Cauchy  relation  is  introduced  in 
order  to  provide  frequency-dependent  expressions  for  the  indices  of  refraction 
and  the  Glads t one-Dale  constant.  Therefore  we  write: 


i) 

a 


1  =  A  +  B  v.2  -t-  C  v4  + . ..... . . 

a  a  a 


K  =  K  +  K  v2  +  Ki  v4+.. ..... 

o  2  4 

At  this  point,  it  is  convenient  and.  instructive  to  introduce 
a  change  of' variable.  Therefore,  let  (Ref.  p): 


where : 


v  +  v 

y  ZZ  Stm  . 


Therefore ; 
When: 


v  - 


V  -r  A  =  v  (  1  -4  ~-  j  ;  dv 
V 


-  dA 


V  -  V 

2 

V  =  V.-,. 


A  =  v2 


_  vz~ 


A  =  v,  -  v  =  -  ysiil. 


We  must  also  rewrite  the  power  series  describing  the  indices  of 
refraction  in  terms  of  the  new  variable.  It  can  oe-  shown  that  (Ref.  5): 


Ai 


-*f  - 

I 


I 


where : 


n  =■  1  +  A  +  B  v2*-  C  v4  + 
a  a  a  a 


■ 1  *  Ai +  A; (4)+  * (4)2  +  Ai(-f)= 


A '  =  A  +  B  “2  +  C  v4  +  D  v6  + 
0  a  a  a  a 


A.  =  2  B  v2  +  4C  v4  +  6  D  V 


A2 


A3 


B  v2  +  6c  V4  +  15  B  v  + 
a  a  a 


_  _A 

4C  v*4  +  20  D  v  + 
a  a 


(A.l) 


Similarly: 


where : 


etc . 


K  =  Kq  +  K2  v2  +  v4 


*  "b +  ^  (-#•) +  a,  (4) +  Aa  (4) 


*t>  *  S0  *  K2  "2  +  h  V4  +  %  -J6  + 


Al' 
a2  = 


2K2  v2  +  4K4  6  K6  T6  +  . . 

K0  v2  +  6K^  "v4  +  15  K6  v6  + 

4K,  T4  +  20  K,  T6  + 

4  o 


?*  (A. 2) 


etc. 

Likewise,  the  term  in  Eq.  43  which  gives  the  contribution  of  the  electrons 
must  be  rewritten: 


L(u  )2 
P  2 

87tvc 


L(w  )2 

P  2 

8ttc~ 


1  + 


f 

!  i 


tvS 

«;•  .f 


r 

.1 


«v: 


2  /  A  /  A  V 

F  C1  -  “  ’  (t> . ) 

4  x3  7  C1  -  4- +  (4)  - . ) 


where:  L^)^  Le^Ng 

X3  ~  dircV*  ~  SttcT^c  m 

o  e 


(A.3) 


A2 


I 


Therefore,  Eq.  43  becomes  finally: 


2 


A  2  +*&- v: 

I  dA  =  f  2  cos' 

Vi  2  .  Vo-  v. 


2  |  Xt  -(p2-pi)(a0  *(4}-~Xl  +  4-) 
-X2  vy(l  ♦  ^  +  Ai  (-f)  +....)(l+  4) 
-x^1’ ?+(£)•  •••)}“ 


Dropping  all  terms  of  second  and  higher  order,  we  obtain: 
iSZJH  A  2  +  .v.gi-v.;.  f  f 

I-  vjv,  r  dA  “  -f—  /_  *  <k»24  v(pp-p.) U  +  4 


Vg-  v1 

L  Vg-Vj  COs2{Xl  7(P2-Pl)[^b  +  (A0  + 

2 

-  X2 "  4  -  x3 7  (x  -  4)}  “  (A-1,) 


This  expression  is  of  the  form: 


J cos2(M  +  Nx)  dx  =  |*  J  cos  2  (M  +  Nx)dx  +  $  J  dx 
=  £jj  sin  2  (M  +  Nx)  +  | 


Thus,  Eq.  A. 4  becomes 
.  2 


I(V2’Vl)  2  7  [x1(p£-n)(A0+  A!)-x2y(i+  A^+A')  +x3]  [sin{2Xi  (p£-Pi)tV 

(Aq+Ai)(^2|^)J  -2X2y  v  +  A£  v  +  (l+A’  +  Ai)(  -  2X3v  +  2 

sin  ■|^2X1(p£-p1)^  AQV  -  (A0+Ax)(^f^)  -2X2y  W-  A'v  -(l+  A^+Al)( 

2X3  ?  -  2X3  (^)j  +  \  (vg-vx)  (A.5) 


(A.5) 


(vg-vx)  =  ^  (vg-VijJ  1  f  aXi(p£-Px)(VA1).X2y(l^  A1>X33 . Mil). 

4  ^  [  X1(p2-p1)(A0+Ai)-X2y(l+A^+  Ai)+X3]  (vg-v*) 

x  cos{  2  [X1(p2-p1)  &0  -  X2y  (1  +  Aq)  -  X^  V  )  j-  (A.6) 


where  Eq.  A.6  is  the  same  as  Eq.  44  in  Sec.  2.2.4. 


As  noted  in  Sec.  2.2.4,  the  phenomenon  of  "rollover"  occurs  because 
the  argument  of  the  sine  contains  factors  (Aq+A^)  and  (l  +  Aq  +  A|),  while  the 


■3 


1  v*?>v'-<»iScc 


i 


>-4*a  ye 


'*'  .3%  V-w*> 


corresponding  factors  in  the  argument  of  the  cosine  are  Aq  and  (l  +  Ai)  .  It 
is  noted,  however,  that  the  effect  is  caused  not  so  much  by  the  fact  that 
Aq  /  Aq  +  A^,  and  (i  +  A^)  j-  (l  *  A^  +  A’ ) ,  but  by  the  fact  that: 


+A 


1 


/ 


1  + 


A0 


i  + 


(A. 7) 


where  it  has  been  assumed,  as  before,  that  the  electron  number  density,  i.e,, 
X^*  is  zero. 

To  show  the  significance  of  these  quantities,  the  group  index  of 
refraction,  ng,  must  be  introduced,  where: 


and 


by  definition. 


cr  V 

0  g 

v  =  =  (^Y1 

g  dk  \du J 

Noting  that; 


;  dw  =  27niv 


Now: 


d  v  =  - 
P 


d  n 

P 


and  at  frequency  v; 


-  n  - 


P 


cv 


(A.8) 


(A. 9) 


A4 


n 


*■  v 


(A. 9a) 


It 


Now,  from  the  Cauchy  relation: 


(A. 7) 


of 


(A. 8) 


(A.9) 


(A. 9a) 


n  =1+A+°T7£+C  v4*. 

PQ  a  a  a 

Ct 


one  obtains  by  differentiation: 

d  nT, 

Pa 


=  il  B7  +  If  C  v 
dV  a  a 


773 


so  that  by  substituting  in  Eq.  A.9a>  one  obtains: 


n6a 


1  +  A  +  B  v2  +  C  v4  + . +  2B  v2+  4  C  v4 

a  a 


a  a  a 

**"“P  r-  «  ““A 
I**  A 


=  1  +  A  +  3B  v“t  5  c  v4+. 
a  a  a 


(A. 10) 


(A. 11) 


Going  back  to  Eq.  A.l,  it  is  seen  that: 


and 


1  +  A '  =  1  +  A  +  B  v2+  C  v4+  ..... . 

0  a  a  a 

1  +  A1  +  A 1  =  1  +  A  +  3B  vz  +  9  C  v4  + 
0  1  a  a  a 


Thus,  it  is  seen  that  to  use  (l  +  A*)  is  to  use  the  phase  index  of  refraction, 
and  to  use  (l  +  A^  +  Al)  is  to  use  the  group  index  of  refraction. 

To  obtain,  the  significance  of  Aq  and  (Aq  +  A-^) ,  one  must  go  back 
to  the  definition  of  the  Gladstone-Dale  constant.  Again  neglecting  the  electrons, 
we  obtain: 

r-  -  Uli 
K  po 

where  pQ  is  the  density  at  which  n  is  measured.  Thus  we  may  vrite: 

....  =  —  (A.  +  3,.  vs  +  C.  v4  + . ) 

Pr,  t  t  z  ’ 


K  =  KQ  +  K2  v2  +  K4  v4  + 


so  that: 


K  =  — 

0  po 

Bt 

K  =  -- 

42  Po 


K,. 


etc. 

and  going  back  to  Eq.  A. 2,  it  is  seen  that: 

A0  =  Ko  +  K2  ‘2  +  \  ‘4  +  . 


+  h  v2  +  Si  -4  + 


A5 


and: 


A0+A, 


V 

is. 


v2+ 


V2+ 


5  K4  “4+ 


ngt  - 1 
Po 


Therefore,  it  is  seen  that  to  use  Aq  is  to  use  the  Gladstone-Dale  constant  based 
on  phase  index  of  refraction,  and  to  use  (Aq  +A,),  is  to  use  the  Gladstone-Dale 
constant  based  on  group  index  of  refraction.  Thus,  substituting  in  Eq.  A, 7,  it 
is  found  that: 


1 


K  B4*  C+ 

7T  +  7T  72  +  7T  "4  +  i-. 

Po..  PQ  .  Po _  ^ 

A.  K 

—  +  3  — 7*4-5  iv4  + . 

P  J  P  P 

o  Ko  Ko 

2B+  v2+  4  C*  v4+  . 

^ _ t _  ^ 

At  +  3Bt  v2+  5  Ct  v4  +  . 


1  +  A  +  B  T2  +  C  -r4  + 
a  a  a 

1  +  A  +  3B  T?2+  5  C  v4  + 
a  a  a 


2  B  ?  +  l*C  v4+ . 

a _ a _ 

1  +  A  +3B  v2  +  5  C  T4  + 
a  a  a 


°r*  A.  +  B,  v2  +  Q.  v4  +  .....  +  2  B.  v2  +  4  C.  v4  + . 

t  t  t  t  t 

/l  +  A+B  v2  +  C  v4  +.,..+  2  B  v2  +  4  C  v4+.... 
'  a  a  a  a  a 

where  we  have  assumed,  for  the  purpose  of  illustration,  that 


A  =  A,  ;  B  =  B,  ;  6  =  C.  ;  etc. 

&  t  £t  "t  Ql  T/ 

Thus  it  can  now  be  seen  why  rollover  occurs;  it  is  because  (l  +  A’)  is  of  the  form 
n  ,  a  number  of  order  one,  and  A_  is  of  the  form  (rip  -  l).  To  both  these  quantities 
ii  added  an  amount  A^(or  A')  which  is  a  number  of  approximately  the  same  order  as 
Aq  or  A' ,  and  A,  forms  a  much  greater  part  of  (A_+  A^than  A!  forms  of  (l+A^+Ap 
so  thatUthe  ratio  Aq/  Aq+A^  is  much  less  than  the  ratio  1  +  Aq/1+Z^  +A£. 

Thus,  in  effect,  rollover  occurs  because  the  Gladstone-Dale  constant  and  the 
index  of  refraction  are  different  quantities,  varying  in  different  ways  with 
frequency.  However,  rollover  is  still  basically  a  dispersive  effect,  since, 
there  was  no  dispersion,  Aj_  and  Aj^  would  be  identically  zero. 

Using  the  data  in  Ref.  6  for  the  index  of  refraction  of  air  and  the 

Gladstone-Dale  constants  for  oxygen,  one  finds  that,  on  fitting  a  Cauchy  relation: 

Kq2  =  (1.843  x  10"4  +  1.773  x  10-35  v2)  -~ 

n_  =  1  +  1.795  x  10-5  +  1.215  x  10"36v2 
*air 

Then,  on  calculating  the  quantities  in  Eq.  A. 6,  one  finds  that,  at  5200$ 

<•  li  m3 

Aq  =  1.902  x  10  i 

A1  =  .1174  x  10*4  m3/kg 


A6 


,  a.  *»»<■***»*»>» 


Aq  =  1.535  x  10"5 
A^  =  .0804  x  ICf5 

Thus  it  is  seen  that  Ai  can  be  neglected  in  relation  to  (l  +  A <j),  since  A|  is  of 
the  order  of  10  ,  while  (l  +  Aq)  is  of  order  unity.  However,  it  is  seen  that 

A-^  is  definitely  not  negligible  with  respect  to  Aq,  and  this  then  leads  to  roll¬ 
over,  as  described  above. 

It  is  instructive,  on  the  basis  of  the  above  data,  to  calculate  the 
fringe  shift  required  for  a  rollover  of  one  fringe  to  occur, i.e.,  for  the  maxi¬ 
mum  of  the  fringe  envelope  to  move  with  respect  to  the  zero-order  fringe  by  an 
amount  equal  to  the  spacing  of  individual  fringes,  so  that  the  maximum  of  the 
fringe  envelope  coincides  with  the  first-order  fringe,  rather  than  the  zero- 
order  fringe. 


Let  y  now  be  the  monochromatic  fringe  shift  at  which  a  rollover  of 
one  occurs.  Since  a  particular  fringe  is  characterized  by  a  particular  value  of 
the  argument  of  the  cosine  in  Eq.  A. 6,  we  may  write,  for  the  zero-order  fringe 
in  particular: 

2  [X1(o2-Pi)  Aq  -  X2y  (1  +  A£)]  v  =  0  (A. 12) 

where  we  have  neglected  electrons.  Substituting  for  X^  and  Xp,  and  putting: 

v  1 


we  obtain  ~(P2-Pi)Aq  -  ^(l  +  Aq)  =  0  (A.13) 

Earlier  it  was  found  that  (l+  Aq)  was  equivalent  to  the  phase  index  of  refraction 
n0o .  Also  from  Sec.  2.2.3,  it  was  found  that  the  spacing  of  individual  fringes 
is  given  by: 

d  - 

2enPa 


Thus,  Eq.  A.13  becomes: 

~~  (P2-Pi)Aq  -  |  =  0  (A.lM 

A  similar  procedure  can  be  carried  out  with  the  argument  of  the 
sine  in  Eq.  A. 6,  which  gives  the  position  of  the  fringe  envelope.  However,  since 
a  rollover  of  one  fringe  has  occurred,  the  position  of  the  maximum  of  the  envelope 
is  at  position  (y  +  d) .  Therefore: 

[X1(o2-Pi)(A0+A1)  -  X2(y  +  d)(l  +  A’)]  (v2-Vl)  =  0 

where  A'  has  been  neglected,  in  line  with  the  order  analysis  carried  out  above. 
Proceeding  as  before,  one  obtains: 

—  (p2-Pi)(Aq+A1)  -  =  0  (A. 15) 
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obtains: 


whence : 


Then,  on  eliminating  =-  (p2~Pi)  between  Eqs.  A.l4  and  A. 15,  one 

(frX*  4  (V  V 


(A.l6) 


Thus,  using  the  data  given  earlier  for  Aq  and  A^,  it  is  found  that,  in  oxygen,  a 
rollover  of  one  fringe  occurs  when  the  monochromatic  fringe  shift  is  about  six¬ 
teen  fringes.  This  agrees  very  well  with  experimentally  observed  rollover  in 
the  oxygen  work. 


If  a  similar  determination  of  Aq  and  A^  is  carried  out  for  argon, 
finds  that  a  rollover  of  one  fringe  occurs  when  the  monochromatic  fringe  shift 
is  about  23  fringes. 


A. 2  White  Light  Interference  with  an  Asymmetric  Intensity  Distribution 

So  far  we  have  considered  only  an  intensity  distribution  which  is 
symmetric  about  V.  We  now  investigate  the  effect  of  an  asymmetric  distribution 
of  intensity,  i.e.,  an  intensity  distribution  with  wavelength  which  is  weighted 
toward  one  end  of  the  wavelength  range,  which  is  the  more  usual  state  of  affairs. 
The  intensity  distribution  in  a  beam  of  white  light,  as  seen  by  the  detector,  is 
a  function  of  many  parameters:  the  spectral  content  of  the  source  light,  the 
transmissions  and  reflectivities  of  the  interferometer  optics,  the  transmissions 
and  reflectivities  of  all  the  other  mirrors  and  lenses  in  the  system,  as  well  as 
the  spectral  response  of  the  detector  itself,  used  to  record  the  interference 
pattern.  Thus,  a  symmetric  intensity  distribution,  such  as  was  considered  above, 
would  be  rare. 


Thus,  referring  to  Fig.  9b,  we  assume  an  intensity  distribution 
with  wavelength  of  the  form: 

**-*?*<*<*  “ 

1=0  for  A  < - ~2cf~  and  A  >  +  * 

O  e id 

In  this  expression,  it  is  noted  that  7  is  a  parameter  which  denotes  the  degree 
of  asymmetry  of  the  intensity  distribution,  and  A  is  the  quantity  introduced 
earlier  which  varies  with  frequency.  Tan” ^7  is,  in  effect,  the  slope  of  the 
top  of  the  intensity  with  frequency  distribution  and  it  is  seen  that  when 
7=0,  this  intensity  distribution  degenerates  to  the  symmetric  distribution 
considered  earlier. 

The  expression  for  the  intensity  distribution  in  the  white-light 
interferogram  then  becomes,  as  in  Eq.  A. 4: 
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The  result  of  the  first  integration  is  identical  to  that  obtained  before  for 
the  symmetric  intensity  distribution,  and  the  result  is  given  in  Eq.  A. 6.  The 
second  integral  is  of  the  form: 

J  x  cos2  (M  +  Nx)  dx  =  -jj-j-  [x  sin  2  (M  +  Nx)  j 

j^cos  2  (M  +  Nx)J  +  "TJ" 


_i_  r 

8m2 


Therefore,  the  second  integral  becomes: 
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(A. 18) 


Thus,  the  complete  expression  for  the  intensity  distribution  in  a  white  light 
interference  pattern  is  the  sum  of  Eqs-.-  A. 6  and  A.l8,  which  is 
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Interpretation  of  Eq.  A.  19  is  very  difficult  to  do  by  inspection,  since  two  terms 
are  involv  ed  which  vary  in  different  ways  with  wavelength.  However,  it  is  noted 
that  Eq.  A.19  reduces  to  Eq.  A. 6  when  7=0. 


M 


A. 3  Numerical  Calculations 


To  illustrate  rollover,  the  effect  of  electrons,  and  the  effect  of 
an  asymmetric  spectral  intensity  distribution,  Eqs.  A. 6  and  A.19  were  programmed 
for  the  University  of  Toronto  IBM  709^  computer,  and  white-light  interference 
patterns  were  calculated  for  several  different  cases  to  illustrate  the  effects 
mentioned  above.  The  conditions,  in  the  final  state  in  the  test-section  were 
assumed  to  be  those  produced  by  a  Mach  12.0  shock  wave  travelling  into  argon  at 
an  initial  pressure  of  21.0  mm  Hg.  This  produced  a  density  change  (Ps-Pi)  of 
about  .175  kg/m3  and  gave  rise  to  about  3*3$  ionization,  i.e.,  sin  electron  number 
density  of  approximately  3«3  x  10^3/m3.  All  the  cases  were  calculated  using 
available  room- temperature  data  for  the  Gladstone-Dale  constants,  and  it  was 
assumed  that  the  Gladstone-Dale  constant  for  the  argon  ion  was  equal  to  that 
for  the  argon  atom.  This  was  done  partly  for  simplicity  (since  these  examples 
were  illustrative  only)  suid  partly  because  the  dispersion  relation  for  the  re¬ 
fractive  index  of  the  ion  (i.e.,  the  variation  of  K^+  with  v)  was  unknown. 


A10 


The  light  from  the  source  which  produced  the  interference  pattern 
was  taken  to  have_a  range  of  wavelengths  from  about  44()0  X  to  5700  A,  with  a 
mean  wavelength,  X,  of  %60  A.  For  the  case  of  y  =  0,  the  intensity  I  (see 
Fig.  9&)  was  evenly  distributed  within  the  above  range;  for  the  case  o?  the 
asymmetric  spectral  intensity  distribution, 7 was  set  ?qual  to  +1.0  (see  Fig.  9b) , 
so  that  the  intensity,  I0,  at  vlt  was  A02/4>  an<i  a-fc  v2  >  *o  was  3A02/4,  so 
that  the  distribution  was  weighted  toward  high  frequencies,  or  short  wavelengths. 

We  consider  first  plots  of  Eq.  A. 6,  i.e., 7  =  0,  so  that  a  spectral 
intensity  distribution  which  is  symmetrical  about  v  is  being  considered. 

Figure  A. la  shows  the  intensity  distribution  of  a  white-light  inter- 
ferogram  with  test  section  and  compensating  chamber  conditions  exactly  equal, 
i.e.,  (P2-P1)  =  0  and  X3  =  0.  It  is  noted  that  the  pattern  is  symmetric  about 
the  position  of  the  zero-order  fringe,  which,  in  this  case,  falls  at  y  =  0,  as 
shown.  Of  particular  note  is  the  phase  inversion  which  occurs  between  the  cen¬ 
tral  fringe  group  and  the  subsidiary  fringe  groups  on  either  side.  This  is 
shown  in  Fig.  A. la,  where  the  order  of  the  fringes  has  been  noted.  It  is  seen 
that  the  position  of  the  4th,  5th,  6th  and  7th  order  fringes  fall  on  minima  in 
the  subsidiary  fringe  groups  rather  than  on  maxima.  Thus,  if  one  wished  to 
measure  fringe  spacing  using  a  fringe  in  the  central  fringe  group  and  a  fringe  in 
the  subsidiary  fringe  group,  one  would  have  to  measure  from  a  maximum  in  the 
central  group  to  a  minimum  in  the  subsidiary  group,  or  vice-versa,  in  order  to 
measure  over  an. integral  number  of  fringes.  It  would  also  appear  that  to  mea¬ 
sure  fringe  spacing  using  two  fringes  one  on  either  side  of  the  same  fringe 
group,  e.g.,  fringes  of  order  +3  and  -3  in  Fig.  A. la,  would  lead  to  small  errors, 
since  the  maxima  of  these  fringes  do  not  fall  exactly  at  the  positions  expected 
on  the  basis  of  the  fringe  spacing,  d,  based  on  v.  This  effect  is  due  to  the 
modulating  influence  of  the  envelope.  A  better  procedure  would  be  to  measure 
from  a  fringe  near  the  centre  of  one  fringe  group  to  a  fringe  near  the  center  of 
an  adjacent  group  taking  into  account  the  phase  inversion  mentioned  above. 

This,  of  course,  assumes  that  the  subsidiary  fringe  groups  are  visible,  which 
may  not  always  be  the  case. 

Figure  A. 2a  shows  the  effect  on  the  white-light  interference 
pattern  of  a  density  change  in  the  test  section,  but  assumes  that  no  electrons 
have  been  produced.  Also  shown  for  reference  is  the  intensity  distribution  on  the 
white-light  interferogram  without  the  density  change,  i.e..  Fig.  A. la.  It  is 
seen  that  the  density  change  causes  both  the  individual  fringes  and  the  fringe 
envelope  to  be  shifted  in  the  direction  of  positive  y,  but  the  envelope  is 
shifted  further  than  the  individual  fringes,  as  shown  by  the  new  positions  of 
the  zero-order  fringe  and  the  maximum  of  the  fringe  envelope.  This  is  the 
"rollover"  effect  discussed  earlier  in  Sec.  2.2.4,  and  is  caused  by  the  fact 
that  the  shift  of  individual  fringes  and  the  shift  of  the  envelope  are  con¬ 
trolled  by  different  indices  of  refraction.  It  is  also  noted  that,  if  the 
density  change  were  large  enough  ,  it  is  possible  for  the  fringe  envelope  to 
be  shifted  sufficiently  far  with  respect  to  the  zero-order  fringe  that  the 
first  order  fringe  might  be  mistaken  for  the  zero-order  fringe.  This  is  the 
danger  that  is  always  present  when  using  white-light  fringes  to  trace  individual 
fringes  through  a  density  discontinuity:  that  the  fringe  which  appears  near 
the  maximum  of  the  fringe  envelope  on  one  side  of  the  discontinuity  may  not  be 
of  the  same  order  as  the  fringe  which  appears  in  the  corresponding  position 
near  the  maximum  of  the  fringe  envelope  on  the  other  side. 


All 


Figure  A. 3a  shows  a  purely  hypothetical  case  in  shock  tube  work: 
here,  electrons  have  been  introduced,  but  no  density  change  has  occurred. 

Again,  Fig.  A. la,  which  gives  the  interference  pattern  before  the  changes  took 
place,  is  shown  for  reference.  It  is  seen  from  the  new  positions  of  the  zero- 
order  fringe  and  of  the  maximum  of  the  fringe  envelope  that  individual  fringes 
have  been  shifted  in  the  direction  of  negative  y  while  the  fringe  envelope 
has  been  shifted  a  similar  amount  toward  positive  y.  It  would  also  appear  from 
this  example  that  the  maximum  of  the  fringe  envelope  coincides  with  the  maxi¬ 
mum  of  one  of  the  fringes;  in  general,  however,  this  will  not  be  true,  and  the 
maximum  of  the  fringe  envelope  would  not  always  coincide  with  the  maximum  of  a 
fringe. 


Finally,  Fig.  A. 4a,  shows  the  effect  of  a  density  change  and  an 
electron  population  appearing  at  the  same  time.  It  is  seen  that  the  effect  of 
rollover,  which  was  illustrated  in  Fig.  A. 2a,  and  the  effect  of  the  electrons, 
which  was  illustrated  in  Fig.  A. 3a,  add  together  to  produce  the  effect  shown 
in  Fig.  A. 4a.  For  example,  it  is  seen  that  the  distance  between  the  positions 
of  the  zero-order  fringe  and  the  maximum  of  the  fringe  envelope  in  Fig.  A. 4a 
is  simply  the  sum  of  the  separations  of  these  two  positions  shown  in  Fig.  A. 2a 
and  Fig.  A. 3a.  Also,  the  shift  of  the  maximum  of  the  fringe  envelope  in  Fig. 

A. 4a  is  the  sum  of  the  shifts  illustrated  in  Fig.  A. 2a  and  Fig.  A. 3a.  The 
same  is  true  for  the  shift  of  the  zero-order  fringe. 

Figure  A. 2a  is  typical  of  the  intensity  distribution  in  an  inter- 
ferogram  of  the  gas  behind -a  shock  wave  in  oxygen,  where  there  is  a  density 
change  and  dissociation,  but  where  the  effect  of  electrons  can  be  neglected. 

Figure  A. 4a  is  typical  of  the  intensity  distribution  in  an  interferogram  of  the 
gas  behind  a  shock  wave  in  argon;  in  this  case  ionization  occurs  giving  rise  to 
an  electron  population,  and  so  producing  the  wide  divergence  between  the  position 
of  the  maximum  of  the  fringe  envelope  and  the  position  of  the  zero-order  fringe. 

In  all  the  above  figures,  the  phase  inversion  between  successive  fringe  groups, 
mentioned  in  connection  .with  Fig.  A. la,  is  apparent. 

Figures  A. la,  A. 2a,  A. 3a,  and  A. 4a  are  plots  of  actual  intensity. 
However,  the  darkening,  or  the  density,  on  a  photographic  plate  is  proportional 
to  the  log  of  the  intensity,  and  this  quantity  is  plotted  in  Figs.  A. lb,  A. 2b, 
and  A. 3h,  and  A. 4b  corresponding  to  the  situations  shown  in  their  companion 
figures.  By  way  of  comparison,  a  typical  densitometer  trace  across  a  white-light 
fringe  pattern  is  shown  in  Fig.  A. 5,  and  it  is  seen  that  the  resemblance  in 
character  between  the  actual  density  distribution  and  the  calculated  distribution 
is  quite  good.  It  is  noted  that  the  minima  are  much  more  sensitive  indicators 
of  the  presence  of  electrons  and  of  rollover  than  are  the  maxima,  so  that  it  is 
often  more  reliable  to  work  with  the  minima  when  attempting  to  estimate  to  what 
extent  these  effects  have  causedthe  position  of  the  maximum  of  the  fringe  envelope 
and  position  of  the  zero-order  fringe  to  separate. 

We  now  turn  to  Eq.  A. 19  where  provision  was  made  for  an  asymmetric 
spectral  intensity  distribution,  as  shown  in  Fig.  9b.  As  mentioned  earlier,  for 
the  calculations  y  was  taken  to  be  +1.0.  Only  three  cases  were  plotted.  Figs. 

A. 6,  A. 7,  and  A. 8,  corresponding  in  all  ways  except  for  the  value  of  y  with  Figs 
A. la,  A. 2a,  and  A. 4a.  It  is  apparent  that  Eq.  A. 19  is  not  a  good  way  of  express¬ 
ing  the  intensity  distribution  in  the  white-light  interferogram,  since,  as  shown 
in  Fig.  A. 6,  the  spacing  of  the  fringes  no  longer  corresponds  to  "F,  but  to  some 
frequency  greater  than  V.  If  y  had  been  chosen  to  be  -1.0,  for  example,  the 
spacing  would  have  corresponded  to  a  frequency  less  than  v.  Plainly  therefore, 

Eq.  A. 19  is  only  approximately  correct,  since  all  the  indices  of  refraction  in 
the  equation  are  expressed  in  terms  of  "v ,  which  is  no  longer  the  frequency 
governing  the  spacing  of  individual  fringes. 
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(Symmetric  spectral  intensity  distribution;  *=  4960A) 


(Symmetric  spectral  intensity  distribution;  *=  496oA) 


(Symmetric  spectral  intensity  distribution;  3  *  4960A 


FIG.A5  -  TYPICAL  DENSITOMETER  TRACE  ACROSS  A  WHITE-LIGHT  FRINGE 
PATTERN  AS  RECORDED  ON  A  PHOTOGRAPHIC  PLATE 


PIG  A6  -  DISTRIBUTION  OF  INTENSITY  IN  A  WHITE- LIGHT  INTERFEROGRAM  WITH  TEST  SECTION  AND 
COMPENSATING  CHAMBER  CONDITIONS  IDENTICAL  -  ASYMMETRIC  SPECTRAL  INTENSITY 
DISTRIBUTION 


APPENDIX  B 


A  Summary  of  the  Experimental  Work  Done  to  Date  (March  1966)  on  the  Determination 

of  the  Gladstone-Pale  Constants  for  Argon 

Some  work  has  been  done  on  the  determination  of  the  Gladstone- 
Dale  constants  for  ionizing  argon,  but  several  problems  have  arisen  preventing 
immediate  completion  of  this  work. 


The  basic  procedure  for  carrying  out  the  determination  of  the 
Gladstone-Dale  constants  for  ionizing  argon  is  given  in  Appendix  B  of  F.ef.  1. 

As  with  the  oxygen  work,  it  was  proposed  to  pass  a  strong  shock  wave  through 
the  gas  and  then  measure  the  index  of  refraction  of  the  resulting  hot  ionizing 
gas  using  the  Mach-Zehnder  interferometer.  Prom  this  determination  of  the 
index  of  refraction,  the  Gladstone-Dale  constants  for  the  argon  atom  and  ion 
can  be  determined,  using  an  anlysis  similar  to  that  outlined  in  Sec. 

There  is,  however,  one  additional  complication  in  the  argon  work: 
electrons  are  present  in  the  hot  gas  as  well  as  ions  and  neutral  atoms,  and  the 
contribution  of  the  electrons  to  the  total  refractivity  of  the  gas  must  be  taken 
into  account.  It  was  proposed  to  do  this  by  making  use  of  simultaneous  white- 
light  and  monochromatic  interferograms .  As  shown  in  Sections  2.2.3,  2.2.4,  and 
Appendix  A,  on  a  given  white-light  interferogram  the  electrons  produce  a  positive 
shift  of  the  white-light  fringe  envelope,  and  a  negative  shift  of  individual 
fringes.  These  shifts  due  to  the  electrons  are  of  equal  magnitude,  and  are  in 
addition  to  say  fringe  shifts  caused  by  a  change  in  density. 

This  can  be  seen  from.Eq.  A. 6.  For  a  given  individual  fringe  on 
a  white-light  interferogram,  the  argument  of  the  cosine  in  Eq.  A. 6  has  a  con¬ 
stant  value.  In  particular,  for  the  zero-order  fringe,  v-  have: 

2  [X1(pi-Pa)A0  -  X2y(l  +  Aq)  -  X^  i  v  =  0  (B.l) 

Since,  as  shown  in  Appendix  A,  Ag  is  the  Gladstone-Dale  constant  corresponding 
to  the  phase  index  of  refraction,  we  substitute: 


V  *0 

Similarly  (-1  +  A')  was  shown  to  be  the  phase  index  of  refraction  of  air;  hence: 

"pa  '  1  +  *0 

Thus,  incorporating  these  changes,  and  substituting  for  X^,  Xg  and  X_,  Eq.  B.l 
becomes:  ** 


£7®V  ,  .  ,  „  ^"Pa  2L7<“P>2 

“ — T—  (P2-P1)  K - -  - - -  =  0 

c  p  c  8ttcv2 


(B.2) 


If  one  now  removes  the  assumption  of  only  one  species  present  in  the  test  section 
(imposed  for  simplicity  in  Sec.  2.2.4),  and  substitutes: 


-v--rr~ 
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*■  ’--.'sv  >■  -  ■*— '•f***' rt»^<**k«sa*  *,.„  , 


,  -  2  -2  7  1 

il7TV  =  t*>  ;  — — •  -  *£-- 


Equation  B.2  becomes 


(“p)l 


f.  .  I  2€npi  L  ^P'2 

?  (Pa-Pi),  KPj_  -  -  y  -  —  -  =  0 

1  J  \>  Jo  2a>2 


From  Sec.  2.2.3,  the  spacing  between  individual  fringes  on  an  interferogram  is 
given  by:  _ 

Ao 


Hence: 


*-  Sip 

d  ■Leim 


2 

=  T  {^[^atorr^1"^  +  KPion  X]  "  Pl  ^atom”  2w2  } 


(B.3) 


where  substitution  has  also  been  made  for  the  Gladstone-Dale  constants;  in  front 
of  the  shock,  only  atoms  are  present,  behind  the  shock,  atoms,  ions  and  electrons 
are  present. 

Similarly,  the  shift  of  the  fringe  envelope  is  given  by  the  argu¬ 
ment  of  the  sine  in  Eq.  A. 6,  vis: 


[X1(p2-Pi)(A0+A1)  -  X2y(l  +  Aq  +  Ai)  +  (vs-vx)  =  0 


(B.4) 


Proceeding,  as  before,  we  substitute  (see  Appendix  A): 

Kg=  Aq  +  Ai 

n8a  =  1  +  ^)  +  ^ 


whence : 


The  approximation  made  here  in  setting  nPa  equal  to  nga  is  in  line  with  the 
note  in  Appendix  A,  where  it  was  pointed  out  that  Aq  and  A£  are  numbers  of 
order  10“4,  and,  since  these  indices  of  refraction  enter  only  through  the 
fringe  spacing,  one  is  justified  in  setting  nPa  equal  to  ng&. 


Thus  proceeding  as  before,  Eq.  B.4  becomes: 


(w  )2 

d  =  Sl2w  =  ”{P£  j^Satora  +  %ion  x  Pl  Kgatom  *  2u^  _} 


(B.5) 


An  expression  similar  to  Eq.  B.3  could  be  written  for  the  case  of 
a  monochromatic  interferogram  also.  For  example,  the  argument  of  cos2  in  Eq.43 


becomes: 


4 

s 


s12m  -  77  {<*  [KPatom  (l-x)  +  KPion  x]  *  p-  Satom 


K) 


2W" 


) 


(B.6) 


where  AQ  in  this  case  is  any  arbitrary  wavelength,  which  may  or  may  not  be 
equal  to  A0  in  Eqs.  B.3  and  B.5. 

_Thus  if  Eqs.  B.3  and  B.5  were  addedp(or  B.5  and  B.6,  provided  A0 
was  equal  to  A0,  i.e.,  w  =  u>),  the  terms  (oj-q) 2/ will  cancel  out,  vis: 


5 12 


m 


s12 


» =  \  i92  - 


'gatom+  ^Patom 


Kl  -  *>  +  (%lon+  KPion>x 


Pl  ^6atom  +  ^Patom^ 


(B.7) 


(Usually,  Spg  will  not  be  visible  on  the  white-light  interferogram,  since  the 
zero  order  fringe  is  too  far  from  the  maximum  of  the  fringe  envelope.  Thus  Si2m 
must  be  obtained  from  a  monochromatic  interferogram,  and  if  Aq  for  the  mono¬ 
chromatic  interferogram  is  not  equal  to  A0  of  the  white-light  interferogram, 
an  error  arises,  because  the  terms  (w  )|/2u)2  will  not  cancel  exactly,  when  Eqs. 
B.5  and  B.6  are  added.)  P 


At  this  stage,  one  must  introduce  an  approximation,  and  set  Kp 
equal  to  Kg.  The  difference  between  these  two  quantities  is  whau  causes  roll¬ 
over,  and,  on  the  basis  of  rollover  measurements,  the  difference  between  Kp  and 
Kg  is  about  5$.  This  reduces  the  number  of.  unknowns  from  four  to  two,  and 
makes  possible  the  use  of  the  analysis  given  in  Sec.  3*1* 

The  preceding  paragraphs  give  the  basis  of  the  method  proposed 
for  determing  the  Gladstone-Dale  constants  for  ionizing  argon.  However  the 
practical  problem  of  determining  the  fringe  shifts  and  ^I2m  rema^ns* 

Once  ionization  begins  behind  a  shock  wave  in  argon,  the  white 
light  fringes  can  no  longer  be  used  to  identify  particular  fringes,  as  was  done 
for  the  oxygen  work,  since*, 'when  ionization  begins  wp2/2ur  is  no  longer  zero 
and  the  shift  of  the  fringe  envelope  on  a  white-light  interferogram  and  the 
shift  of  individual  fringes  on  a  monochromatic  interferogram  are  no  longer  re¬ 
lated  as  they  were  for-.the  oxygen  work.  Hence,  if  the  white-light  fringes  are 
to  be  used  to  trace  particular  individual  fringes  through  the  shock  front,  this 
must  be  done  where  there  is  no  ionization,  i.e.,  immediately  behind  the  shock 
wave.  Once  the  individual  fringes  are  traced  through  the  shock  front,  they 
can  be  followed  visually  to  the  equilibrium  region,,  and  Sp2m  determined. 

Thus,  it  is  seen  that,  if  the  region  immediately  behind  the  shock 
wave  is  obscured  for  some  reason,  individual  fringes  behind  the  shock  wave 
cannot  be  correlated  with  those  in  front  of  the  shock  wave,  and  Si2m  cannot  be 
determined.  This  was  the  first  problem. encountered,  and  is  illustrated  by  the 
intereferograms  shown  in  Fig.  B.l.  These  interferograms  were  taken  with  a 
magnesium  spark  source  in  the  interferometer,  and  the  blurring  of  the  frozen 
(no  ionization)  region  behind  the  shock  wave  is  due  to  the  poor  time  resolution 
(2  psec.)  available  with  this  source.  In  addition,  the  white-light  interfero¬ 
grams  were  badly  fogged  by  the  radiation  of  the  hot  gas  in  the  test  section. 
(Fogging  of  the  monochromatic  interferograms  was  prevented  because  of  an  inter- 
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ferej.ce  filter  between  the  photographic  plate  and  the  test  section  which 
filtered  out  most  of  the  radiation). 

The  expected  fringe  shifts  are  shown  in  Table  El,  along  with  the 
measured  fringe  shifts.  For  the  interferograms  shown  in  Fig.  B.l,  the  frozen 
fringe  shift  and  the  monochromatic  fringe  shift  could  not  be  determined  for  the 
reasons  given  above,  but  the  shift  of  the  white-light  fringe  envelope  appears 
to  be  approximately  correct:  approximately  correct  inasmuch  as  it  is  impossible 
on  these  interferograms  to  determine  exactly  where  the  maximum  of  the  white- 
light  fringe  envelope  falls.  However,  there  appears  to  be  a  reasonable  relaxa¬ 
tion  time  behind  the  shock  wave,  amounting  to  3-4  microseconds.  Thus,  while 
no  accurate  measurements  could  be  made,  the  rough  measurement  of  Sx2y  indicated 
that  conditions  behind  the  shock  wave  were  fairly  close  to  those  predicted  by 
theory.  At  this  point,  development  of  the  exploding  wire  light  source  was  put 
in  hand  so  that  a  Kerr  Cell  could  be  used  to  improve  the  time  resolution  from 
two  microseconds  to  200  nanoseconds,  allowing  accurate  quantitative  information 
to  be  obtained. 


| 

s 
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Some  time  later,  after  development  of  the  exploding  wire  light 
source  was  completed,  a  return  was  made  to  the  argon  work,  and  the  interfero¬ 
grams  from  one  of  the  first  runs  are  shown  in  Fig.  B.2.  It  is  seen  that  the 
required  time  resolution  was  achieved.  Approximate  fringe  shifts  were  deter¬ 
mined  and  are  compared  with  the  predicted  fringe  shifts  in  Table  B.l.  The 
frozen  fringe  shift  was  determined  by  putting  a  piece  of  paper  over  the  flow 
side  of  the  interferogram  leaving  only  the  frozen  region  behind  the  shock  wave 
exposed  to  view.  The  white-light  fringe  pattern  in  the  exposed  frozen  region 
could  then  be  correlated  with  the  white-light  fringe  pattern  in  front  of  the 
shock  wave,  and  particular  fringes  could  then  be  traced  through  the  shock  front, 
as  was  done  for  the  oxygen  work.  The  equilibrium  shift  of  the  white- light 
fringe  envelope  Si2w  was  determined  to  the  nearest  fringe  only,  introducing  an 
error  of  0.5  fringes  at  most.  Si2n  was  determined  by  transferring  the  frozen 
fringe  shift  from  the  white  light  interferogram  to  the  monochromatic  interfero¬ 
gram,  as  was  done  in  the  oxygen  work;  the  individual  fringe  could  then  be  traced 
to’  the  equilibrium  region,  where  Si2m  was  determined. 

These  measurements  showed  that  the  frozen  fringe  shift  was  being 
achieved  roughly  as  predicted,  tout  that  wide  discrepancies  existed  in  S^2m  and- 
Si2w*  These  discrepancies  indicated  that  full  ionization,  as  predicted  by 
normal  shock  calculations,  was  not  being  achieved.  It  was  also  noted  that  the 
relaxation  zone  shown  in  Fig.  B.2  was  much  shorter  than  that  shown  in  Fig.  B.l. 

Some  diagnostic  runs  were  then  taken  in  an* attempt  to  clarify 
the  situation.  Since  the  interferograms  shown  in  Fig.  B.2  were  among  the 
first  taken  with  the  exploding  wire  source  and  Kerr  Cell,  this  run  was  dupli¬ 
cated,  but  using  the  old  magnesium  spark  source  in  order  to  be  certain  that 
the  Kerr  Cell  and/or  exploding  wire  sourc  ■  were  not  causing  odd  effects.  The 
resulting  interferograms  are  shown  in  Fig.  B.3*  The  fogging  due  to  the  test- 
section  radiation  is  apparent,  as  in  Fig.  B.l.  However,  although,  although  the 
frozen  and  monochromatic  fringe  shifts  could  not  be  determined,  the  equili¬ 
brium  shift  of  the  white-light  fringe  envelope* was  as  in  Fig.  B.2  (see  Table 
Dl).  In  addition  Fig.  B.3  indicates  little  or  no  relaxation  zone,  compared  to 
Fig.  B.l,  indicating  that  the  small  relaxation  zone  visible  in  Fig.  B.2  (where 
the  time  resolution  is  better)  is  real.  Thus,  the  run  indicated  that  the 
optical  arrangements  were  not  causing  trouble. 
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At  first  glance  then,  on  the  basis  of  Figs.  B.2  and  B.3>  the  most 
likely  explanation  of  the  anomalies  noted  in  Fig.  B.2  was  an  impurity  in  the  gas 
probably  oxygen  or  nitrogen.  The  dissociation  energy  of  nitrogen  or  oxygen  is 
much  lower  than  the  ionization  energy  of  argon,  so  that  practically  all  of  the 
diatomic  impurity  is  dissociated  before  the  argon  begins  to  ionize.  However, 
only  10$  ionization  was  expected,  so  that,  if  even  only  a  few  percent  of  im¬ 
purity  were  present,  most  of  the  energy  of  the  gas  would  go  into  dissociating 
the  impurity,  and  very  little  into  ionizing  the  argon,  thus  reducing  the  ioniza¬ 
tion  level  below  that  predicted. 

Thus  the  run  illustrated  in  Fig.  B.2  was  duplicated  again,  but 
using  special  high-purity  reagent  grade  argon,  and  the  results  of  this  run  are 
shown  in  Fig.  B.4.  Approximate  measurements  of  fringe  shift  were  made  as  before 
and  these  are  compared  with  theory  in  Table  B.l.  It  is  seen  that  the  results 
are  nearly  identical  with  the  run  shown  in  Fig.  B.2  as  regards  fringe  shifts, 
and  also  as  regards  relaxation  time.  The  conclusion  was  then  reached  that  the 
anomalies  were  not  caused  by  gas  impurity. 


TABLE  B.l 


THEORETICAL  &  MEASURED  FRINGE  SHIFTS  FOR  ARGON  RUNS 
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FIG.  B.l 

Ms  »  13-95, 

=  9.0  mm  Hg,  x  =  8.25$, 

P2/ Pi  =6.18 

Theoretical 

Measured 

s12f 

1.78 

- 

S12w 

6.47 

• 

6.0 

S12m 

-  .50 

- 

FIG.  B.2 

Ms  =  14.29, 

Pj_  =  9.8  mm  Hg,  x  =  8.95$, 

Ps/Px  =6.32 

Theoretical 

mr 

Measured 

s12f 

1.91 

3.4 

S12w 

7.54 

5.5 

S12m 

-  .88 

+2.6 

FIG.  B.3 

Mfi  =  14.18, 

p^  =  9.8  mm  Hg,  x  =  8.65^, 

Ps/Px  ■-  6.26 

Theoretical 

Measured 

s12f 

1.92 

- 

Sl2w 

7.33 

5.5 

S12m 

-  .74 

- 

FIG.  B.4 

Ms  =  l4.l8, 

p^  =  9.8  mm  Hg,  x  =  8.65^, 

P2/P1  =  6.26 

Theoretical 

Measured 

* 

s12f 

1.92 

3.6 

S12w 

7.33 

5.5 

S12m 

-  .74 

+3.0 
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FIG.  Bl  -  WHITE-LIGHT  AND  MONOCHROMATIC  INTERFEROGRAMS  of  a  shock  wave 
in  argon,  typical  of  the  early  interferograms  taken  with  the  magnesium  spark 
light  srurce.  Note  the  blurring  *of  the  shock  front  due  to  poor  time  resolution; 
also  fogging  of  white-light  interferogram  (Mg  =  13. 95  ,  =  9. 0  mm  Hg) 

(Shock  moving  right  to  left;  height  of  field  of  view  is  approximately  7”) 


if 


I 


i 


! 

I 

i 


i 


> 

f 

i 

t 

\ 


t 

! 

i 


fe 


:1 


FIG.  B2  -  WHITE-LIGHT  AND  MONOCHROMATIC  INTERFEROGRAMS  of  a  shock  wave 
in  argon  taken  with  exploding  wire  light  source  and  Kerr  Cell.  Note  improved 
time  resolution  and  lack  of  fogging  (Mg  =  14.29  .  pj  *  9.8  mm  Hg) 

(Shock  moving  right  to  left;  height  of  field  of  view  is  approximately  7") 
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APPENDIX  C 

An  Exploding  Wire  Light  Source  for  the  UTIAS  9- in.  Mach-Zehnder  Interferometer' 
C.l  Requirement 

In  the  course  of  the  experimental  work  to  determine  the  Gladstone-' 
Dale  constants,  a  need  arose  for  improved  time  resolution  in  the  interferograms 
of  the  shock  wave.  Previously*  the  light  source  used  with  the  interferometer 
had  been  a  spark  using  magnesium  electrodes,  with  energy  supplied  by  a 
0.2  uf  capacitor  charged  to  10  kv.  A  typical  time-resolved  intensity  trace  of 
the  emitted  light,  take»  with  an  EG  &  G*  SD-100  photodiode,  is  shown  in  Fig.  1: 
it  is  seen  that  the  duration  of  the  light  output  is  of  the  order  of  two  micro¬ 
seconds  . 


This  was  sufficiently  brief  for  the  trial's  "using  oxygen  as'  the  ' 
test  gas,  even  though  the  shock  wave  was  blurred  considerably,  since  it  was 
necessary  only  that  the  equilibrium  regions  be  clearly  defined.  For  the  argOn 
work,  however,  this  was  no  longer  acceptable.  Here  it  was  necessary  that  the 
Rankine-Hugoniot  region  immediately  behind  the  shock  wave,  where  there  was  no 
ionization,  be  visible,  (see  Appendix  B)  and  thus  time  resolution  of  better  than 
one  microsecond  became  necessary. 

The  increased  time  resolution  could  be  obtained  in  either  of' two 
ways:  by  decreasing  the  duration  of  the  light  output  from  the  spark,  or  by  in¬ 
stalling  in  the  system  a  very  fast  shutter,  in  either  case,  cutting  down  the 
duration  of  the  light  output  as  seen  by  the  camera.  Both  of  these  methods  were, 
investigated.  It  was  found  however  that,  to  shorten  the  duration  of  the  spark, 
it  was  necessary  to  reduce  the  energy  of  the  discharge,  and  no  significant 
decrease  in  duration  occurred  until  the  energy  of  the  discharge  was  so  low  .that 
the  resultant  light  output  was  too  weak  for  photographic  recording..  Various 
kinds  of  capacitors  in  various  configurations  were  investigated*  but  without 
significant  improvement. 

Thus,  only  the  alternative  remained.  However,  no  mechanical 
shutter  could  be  made  fast  enough  to  satisfy  the  requirements  for  time;  resolution 
and  still  be  capable  of  precise  synchronization  with  other  events.  This,  then, 
necessitated  the  use  of  an  electro-optical  shutter,  e.g.,  a  Kerr  Cell,  which  has 
no  moving  parts  and  can  be  synchronized  precisely,  Unfortunately,  a  Kerr  Cell 
has;  inherently,  very  low  transmission**even  when  open,  and  it  was  found,  when 
one  was  used  in  conjunction  with  the  magnesium  spark,  that  the  intensity  of  the 
light  output  was  too  low  to  overcome  the  poor  transmission.  Efforts  to  increase 
the  intensity  of  the  light  from  the  source  by  increasing  the  energy  of  the  dis¬ 
charge  produced  very  little  improvement  and  caused  excessively  fast  erosion  of 
the  magnesium  electrodes. 


*  Edgerton,  Germeshausen  &  Grier,  Inc., 


**  The  theoretical  maximum  transmission  of  a  beam  of  non- polarized  light  is  only 
5036;  only  the  use  of  expensive  crystal  polarizers  will  allow  the  transmission  to 
approach  this  figure. 


However ,  with  the  Kerr  Cell,  the  duration  of  the  light  output 
from  the  source  became  of  secondary  importance  end  brevity  could  be  sacrificed 
for  increased  peak  intensity .  Thus  a  different^ type  of  light  source  could  be 
used;;insplacerof xthe^epwrkv 

The  requirement  was  for  a  source  of  small  diameter*  since  the 
interferometer  required  at  least  k  pseudo-point  source,  confined  with  high  peak 
brightness,  i.e.,  light  output  per  unit  urea.  For  this  reason,  flash  tubes  and 
other  types  of  lamps  were  M4ect^»  thebe  had  plenty  of  power,  but  their 
luminous  area  was  quite  large,  leading  to  low  brightness.  Other  light  sources 
involving  explosives  were  :reject(»d  because  of  possible  damage  to  the  interfero¬ 
meter  situated  close  by; 

This  left  the  possibility  of  using  the  light  output,  obtained 
when  a  small  piece  of  wire  was  vaporized  by  an  electric  discharge.  Our  coh- 
.  tacts  with  other  laboratories  (Ref.  16)  had  indicated  that  the  exploding  wire 

solution*.  having  been  used  with  considerable  success. 

TOere  remaihed^ati  this  stage,  one  further,  possibility:  that 
of  Using  a  laser  as  the- light  Source.  This,  however,  would^have:  involved  con¬ 
siderable  expense  and  delay  besides  being  difficult  to  adapt  to  the  present 
installation.  Thus  the  idea  of  using  a  laser  was  rejected,  at  least  for  the 
time  being,.  Since  the  exploding  wire  source  appeared  to  be  much  cheaper,  easier 
to  construct,  and  more  ejq*dient  in  terms  of  the  time  required  to  get  it  into 
J .operatiioni,  s "  "" 

CV2  Factors  Controlling  Design  *  .* 

In  the  des^  «^  the  ej^l<>ding  wire  light  source,  a  number  of 
criteria  had  tbbf  kepi  dn  miiSf;  V  '  ” 

1)  The  size  and  composition  of  the  wire  to  be  vaporized.  Reference 

l6  indicated  that  tungsten  wiire,  .001-in.  diameter  and.  approximately  1/4"  long, 
would  be  suitable.  This  Was  taken  as  a  starting:  point,  and  the  #est  of  the  unit 
was  designed  around  this.  i^proation^  " 

2)  Availability  of  components.  Some  kind  of  switch  was  required  to 
trigger  the  discharge  at  a  precise  instant,  add  there  was  available  an  EG  &  6 
Model  GP-12A  Triggered  Spark  Gap  Which  appeared  suitable  for  this  purpose:  the 
unit  was  sufficiently  rugged  to  withstand  the  massive  electrical  discharges 
involved,  and  was  capable  of  precise  triggering.-'  A  large  capacitor  was  also 
required,  and,  in  line  with  the  suggestion  of  Ref;  16,  there  was  available  a 
BICC*  8.5l*f,  20  kv  Energy  Storage  Capacitor,  capable  of  storing  up  to  1600 
joules.  Thus  the  physical  dimensions  and  characteristics  of  these  two  items 
had  to.  be  kept  in  mind  during  des.ign.^ 

3)  Installation  requirements.  In  order  to  provide  reasonable  inter¬ 
changeability  between  the  exploding  wire  source  and  the  magnesium  spark  source, 
as  much  as  possible  of  the  equipment  associated  with  the  existing  light  source 
installation  had  to  be  used.  This  meant  using  the'  existing  light  source  mount¬ 
ing,  which  put  rather  severe  limitations  oh  the  physical  layout  of  the  compon¬ 
ents,  and  dictated,  among  other  things,  that  the  capacitor  had  to  sit  on  the 

_  -  -  >  3  ,v  >  * 

f  British  Insulated  Callender's  Cables  Limited. 
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floor  outside  the  light  source  box,  as  shown  in  Fig.  C.3»  and  be  connected  to 
it:  ;by  'd'ffrexibre  cable .  This  requirement  also  demanded  the  use  of  the  exist* 
ing  :high-voltage  power  supply,  and,  as  much  as  possible,  the  existing  trigger 
spark:  generator  used  to  fire  the  triggered  spark  gap. 

U)  High  voltage  insulation.  Given  the  physical  dimensions  and  other 

characteristics  of  the  components  and  the  installation  requirements,  due  regard 
had  to  be  given  to  high  voltage  insulation.  This  becomes  especially  important 
when  one  considers  the  severe  and  extensive  damage  that  could  result  to  the 
interferometer  and  its  valuable  optics,  situated  close  by,  if  a  breakdown  in 
the- insulation  occurred  with  the  capacitor  fully  charged  to  1600  joules. 

5)  Finally,  given  all  the  above  requirements, it  was  desirable  to 
keep > inductance  ^uid.  resistance  in  the  electrical  circuit  as  low  as  possible. 

*  Low  inductance-  Was  desirable  in  order  to  get  a  fast  risetime  for  the  energy 
discharge;,  this  would  help  decrease  the  duration  and  increase  the  peak  in¬ 
tensity  pf  the  light  output.  This  requirement  dictated,  for  examble,  the  use 
of.  coaxial  vcable  to  connect  the  capacitor  to  the  source.  Low  resistance  was 
desirable  in  order  to  keep  ohmic  losses  to  a  minimum  so  that  as  much  energy 
as -possible  would  go  into  vaporizing  the  wire.  Even  a  few  ohms  are  significant 
here,  as:  peak  currents  are  of  the  order  of  lCr  amperes. 

6)  In  addition,  it  was  desired  that  use  of  the  exploding  wire  source 
-  should  be  as  convenient  as  possible,  with  low  "turn-around"  time  between 

successive  shots ; .  Given  the  installation  requirements,  this  demanded  that  the 
light  source  unit  proper  be  easily  removable  in  order  to  replace  the  wire. 


Before  proceeding  with  an  actual  design  of  an  exploding  wire 
sourcer  for  use  :withvthe  :iht^i%erometer,  a  preliminary  model  was  built  in  order 
to  -.test  the  concept,,  and  get : some  idea  of .  the  performance  that  might  be  expect¬ 
ed^  This,  model;  was  quite  ..rudimentary,  designed  only  to  deliver  the  charge  of 
the  capacitor  through:  the  triggered  spark  gap  to  the  requisite  piece  of  tungs- 
ten  wire.  Triggering.  te.chniques  for  the  spark  .gap  were  also  developed  using 
this  model, .-with'  the,  result  that  an  EG  &  G  Model  TR-60  Pulse  Transformer  was 
*  .  siRjstituted  for  the  existing- transformer  in  the  trigger  spark  generator. 

Polarity  requirements  for  the  trigger  spark  were  also  determined. 

The  light,  output  from  this  arrangement  was  monitored,  as  was  the 
magnesium  spark,-  by-  an  EG&G  Model  SD?100  Photodiode,  which  had  rise  and  fall 
timesofHthe^order-of  fiyeand  fifteen  nanoseconds  respectively,  so  that  the 
:Intehsity;profilewas  faithfully  reproduced  without  complications  caused  by 
•the,  sensor j  Intensity  measurements  with  this  arrangement  indicated  that  an 
improvement  in  peak  intensity  of  at  least  3  -  4  times  over  that  of  the  mag¬ 
nesium  spark  couid  be  expected.  Thus,  it  was  considered  that,  along  with  other 
improvements  which  could  te  made  in  the  system  as  a  whole,  an  exploding  wire 
.  source  pluus  aKerr  CeLL  Shutter  would  be  a  practical  way  of  obtaining  the  requir¬ 
ed  time  resolution . 

=.  ,  •  Design  was  then  put  in  hand  for  a  final  model  to  be  used  with 

the.  interferometer and  drawings  and  photographs  are  shown  in  Figs.-  C.3  -  C.6. 
The>;bas^ic,i\^'tjT-as~shoTO/ih7Fig;C.h,  is  built  of  Plexiglas  and  is  almost  totally 
ehcl08ed  fdr  high  voltage  insulation.  This  unit  sits  in  a  Plexiglas  cradle, 
visible  in  Fig..  3,.  which  is  boited  to  the  light  -source  box  in  place  of  the 
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original  magnesium  spark  unit,  so  that  the  two  light  sources  are  interchangeable. 
As  mentioned  before,  the  capacitor  had  to  be  placed  outside  the  light  source  on 
the  floor  and  connected  to  it  by  cable.  For  this  purpose,  a  short  length  of 
RG-17A/U  coaxial  cable,  rated  at  11,000  volts  at  radio  frequencies  (higher  at 
lower  frequencies),  was  obtained;  this  cable  was  reasonably  flexible,  yet  had 
sufficiently  big  conductors  that  resistive  losses  were  minimized.  The  termina¬ 
tion  of  the  cable  at  the  capacitor  was  built  of  brass,  and  its  design  was  dic¬ 
tated  by  the  form  of  the  connections  on  top  of  the  capacitor.  At  the  other  end, 
in  order  that  the  light  source  unit  should  be  quickly  removable,  the  use  of 
some  kind  of  screw-on  connector  was  considered  desireable;  however,  a  connector 
big  enough  to  withstand  the  currents  involved,  and  with  sufficient  stand-off 
distance  to  hold  off  20,000  volts  was  not  readily  available.  Eventually,  a 
large  Type  LC  radio- frequency  connector  was  obtained  (Plug:  UG-15VU> 

Jack:  UG-352/U)  rated  at  only  5000  volts,  but  sufficiently  robust  otherwise 
for  this  application.  A  new  teflon  insert  was  then  made  for  the  jack,  and 
the  assembly  of  the  plug  was  modified  slightly,  in  order  to  increase  the  stand¬ 
off  distance.  Details  are  shown  in  Fig.  C.6. 
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The  electrical  circuit  is  shown  in  Fig.  C.7.  It  is  noted  that 
the  high-voltage  connection  from  the  power  supply  is  made  to  the  top  electrode 
of  the  spark  gap,  that  a  twenty  megohm,  forty  watt  resistor  is  installed  in  the 
charging  line  (to  protect  the  power  supply  from  the  discharge),  and  that  the 
light  source  is  grounded  on  the  ground  side  of  the  capacitor  rather  than  between 
the  spark  gap  and  the  tungsten  wire.  This  last  was  to  ensure  that  all  the 
energy  stored  in  the  capacitor  went  through  the  tungsten  wire  and  also  so  that 
the  discharge  would  not  momentarily  raise  the  ground  point  above  true  ground 
potential,  affecting  all  of  the  instrumentation  connected  to  it.  This  effect 
will  always  occur  to  a  certain  extent,  but  the  resistance  between  the  ground 
point  and  the  capacitor  is  less  in  this  configuration  than  if  the  ground 
connection  had  been  between  spark  gap  and  tungsten  wire.  The  advantage  of  the 
latter  configuration  would  be  that  the  energy  of  the  trigger  spark  would  not 
pass  through  the  wire  to  be  exploded;  as  it  stands,  the  trigger  spark  energy 
does  pass  through  the  wire,  and  this  sets  a  lower  limit  to  the  diameter  of 
the  wire:  if  the  wire  is  too  fine,  the  energy  of  the  trigger  spark  will 
vaporize  it  before  the  main  discharge  begins. 

The  control  and  monitoring  circuitry  is  also  shown  in  Fig.  C.7, 
and  photographs  of  the  entire  installation  and  the  Kerr  Cell  installation  are 
shown  in  Fig.  C.8.  The  entire  arrangement  is  triggered  by  a  +150  volt  signal 
from  a  pulse  generator  which  is,  in  turn,  fired  by  a  shock  detector  located  in 
the  wall  of  the  shock  tube  a  few  feet  upstream  of  the  test  section.  This  pri¬ 
mary  150  volt  trigger  is  used  to  fire  the  first  delay  generator,  a  Racal  digital 
delay  generator,  with  the  delay  set  to  prevent  any  further  operation  until  the 
shock  wave  or  other  phenomenon  is  in  a  position  to  be  photographed.  The  first 
delay  generator  then  puts  out  a  +20  volt  signal  which  fires  another  pulse  genera¬ 
tor,  which  puts  out  a  secondary  +150  volt  trigger  signal.  This  secondary 
trigger  is  then  attenuated  and  used  to  fire  the  trigger  spark  generator  which 
provides  the  pulse  for  initiation  of  the  main  discharge.  The  secondary  trigger 
is  also  used  to  fire  a  second  delay  generator,  an  Electro-Optical  Instruments 
Model  TG-100  Delay  Generator,  which  has  a  continuously  variable  time  delay 
extending  down  to  about  100  nanoseconds.  This  second  delay  generator  is  used 
to  delay  the  trigger  pulse  to  the  Kerr  Cell  until  the  light  output  from  the 
exploding  wire  has  reached  peak  intensity,  approximately  two  micro-seconds 
after  initiation  of  the  discharge.  The  output  from  the  first  delay  generator 
is  also  used  to  trigger  an  oscilloscope  which  records  the  signal  from  an 
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*  SD-100  photodiode  mounted  inside  the  interferometer  (in  a  position  to  watch  the 
light  output  from  the  exploding  wire)  and  also  records  a  sample  of  the 
electrical  pulse  which  operates  the  Kerr  Cell  (See  Fig.  C.ll.)  These  two  are 
recorded  simultaneously,  and  thus  the  synchronization  between  the  Kerr  Cell 
and  the  discharge  can  be  checked. 

C.4  Performance 


Figure  C.2  shows  a  typical  intensity  profile  obtained  using  the 
SD-100  photodiode.  Since  the  traces  shown  in  Figs.  C.2  and  C.l  were  obtained 
under  identical  conditions,  it  is  seen  that  an  improvement  in  peak  intensity 
of  approximately  six  times  over  that  of  the  magnesium  spark  was  achieved, 
assuming  that  the  response  of  the  photodiode  was  linear  with  incident  intensity, 
something  which  was  open  to  doubt:  a  similar  trial  using  a  different  SD-100 
photodiode  showed  only  3-1/2  -  4  times  improvement  over  the  magnesium  spark.  In 
addition  it  was  noted  that  the  total  duration  of  the  light  had  increased  to 
approximately  twenty  microseconds.  Intensity  traces  were  also  obtained  using 
three  different  interference  filters  as  shown  in  Fig.  C.$:  it  is  seen  that  the 
profiles  tended  to  vary  somewhat  with  wavelength,  something  which  was  not  un¬ 
expected.  Again,  the  peak  intensity  appeared  to  vary,  but  this  was  most  likely 
due  to  the  varying  response  with  wavelength  of  the  photodiode,  as  well  as  being 
due  to  the  varying  transmission  of  the  filters. 

A  spectral  analysis  of  the  light  output  was  also  attempted,  but, 
apart  from  providing  some  indication  of  the  combined  characteristics  of  the 
light  source,  interferometer  optics,  Kerr  Cell  (if  included),  spectrograph  and 
film,  little  was  accomplished.  Sample  spectra  are  shown  in  Fig.  C.10.  The 
main  thing  that  was  noted  from  these  was  the  cutoff  of  the  Kerr  Cell,  at  around 
440QA.  It  proved  impossible  to  correlate  the  line's  of  these  spectra  with  the 
lines  of  a  tungsten  spectrum  with  any  degree  of  certainty.  This  was  not  too 
surprising  as,  in  addition  to  tungsten  lines,  there  were  undoubtedly  lines  pre¬ 
sent  from  all  the  brass  and  plastic  parts  in  the  source  adjacent  to  the  explod¬ 
ing  wire,  as  well  as  from  the  air.  The  relative  intensities  of  the  lines  were 
also  distorted,  since  the  dispersion  pf  the  spectrograph  was  not  constant  with 
wavelength,  and  since  the  spectra  were  taken  almost  at  the  limit  of  sensitivity 
of  the  film,  with  a  very  wide  slit  setting  on  the  spectrograph.  In  addition, 
the  lines  were  very  likely  broadened  and.  possibly  shifted  from  their  correct 
positions  by  the  intense  electric  and  magnetic  fields  present  in  the  region  of 
the  wire. 


One  positive  result  of  this  analysis  was  to  emphasi that, 
while  it  was  now  possible  to  use  photographic  film  for  interferometry  using 
the  exploding  wire  source  plus  the  Kerr  Cell,  the  performance  was  still  mar¬ 
ginal.  Therefore,  some  experiments* were  carried  out  to  find  the  combination 
of  photographic  plate  plus  developer  that  would  produce  the  best  results. 

After  trying  several  combinations,  it  was  judged  that  Ilford  HPS  plates  develop¬ 
ed  for  20-30  minutes  in  May  and  Baker  "Promicrol"  developer  proc..?<.d  the  best 
results*:  as  good  contrast  as  possible,  combined  with  low  fog  level  and  rea¬ 
sonably  fine  grain.  Even  then,  a  considerable  reduction  in  the  size  of  the 
image  was  necessary  in  order  to  produce  useable  interferograms.  However,  the 
reduction  in  size  (from  approximately  4"  x  3-1/4"  to  1-3/4"  x  1-1/2")  did  not 
prove  too  seriops  in  terms  of  the  accuracy  with  which  measurements  of  fringe 
shift  could  be  made  on  the  plates.  However,  even  with  the  reduction  in  size, 

*  Latest  trials  have  indicated  that  Kodak  2475  emulsion,  donated  by  Eastman 
Kodak,  Rochester,  may  offer  significant  improvement. 
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this  was  a  great  improvement  over  the  magnesium  spark  plus  Kerr  Cell  combination, 
where  it  proved  to  be  impossible  to  record  on  film  under  any  circumstances. 

In  general,  it  was  found  that  a  number  of  things  had  to  be 
watched  carefully.  One  of  the  main  sources  of  trouble  was  feedback  from  the 
discharge  through  the  trigger  spark  generator  to  the  oscilloscope,  and  particu¬ 
larly  to  the  second  delay  generator.  This  feedback  was  of  considerable  magni¬ 
tude  and  would  cause  the  second  delay  generator,  and  thus  the  Kerr  Cell,  to 
trigger  prematurely.  It  was  this  feedback,  plus  the  fact  that  the  second  delay 
generator  required  a  pulse  of  about  100  volts  to  trigger  consistently,  that 
forced  the  use  of  the  control  system  shown  in  Fig.  C.7,  where  the  output  of 
the  first  delay  generator  was  used  to  fire  a  pulse  generator  which  provided  a 
secondary  trigger  pulse.  This  secondary  trigger  had  to  be  attenuated  before 
being  used  to  fire  the  trigger  spark  generator,  but  the  same  attenuator  also 
cut  down  any  feedback  coming  from  the  discharge  to  the  point  where  it  would 
not  affect  the  remaining  instrumentation.  Originally,  it  had  been  hoped  to 
fire  the  second  delay  generator  from  the  discharge  itself,  an  arrangement  which 
would  have  nullified  any  jitter  in  the  triggering  of  the  main  discharge.  As  it 
stands,  in  the  configuration  shown  in  Fig.  C.7,  there  is  a  small  amount  of 
jitter,  about  0.3  microseconds,  in  the  triggering  of  the  discharge.  However, 
as  seen  in  Figs.  C.2  and  C.9,  the  peak  of  the  intensity  curve  is  broad  enough 
that  this  amount  of  jitter  can  be  tolerated.  In  addition,  the  synchronization 
between  the  Kerr  Cell  and  discharge  was  monitored  with  every  shot,  so  that  minor 
drifts  in  the  triggering  levels  of  the  various  components  could  be  spotted  and 
compensated  for . 

Figure  C.ll  shows  a  typical  oscilloscope  trace  taken  to  check 
.  synchronization.  The  discharge  is  monitored  (as  mentioned  earlier)  by  a  photo¬ 
diode  placed  inside  the  interferometer  and  considerably  closer  to  the  source 
than  was  the  case  in  Figs.  C.2  and  C.9.  Hence,  the  electrical  pickup  from  the 
discharge  is  much  Horse.  The  traces  shown  in  Figs.  C.2  and  C.9  were  taken  with 
the  photodiode  placed  in  the  position  normally  occupied  by  the  camera,  something 
which  is,  of  course,  not  possible  when  taking  interf erograms .  The  upper  trace, 
as  shown  in  Fig..  C.ll,  records  a  sample  cf  the  35,000  volts  pulse  which  operates 
the  Kerr  Cell.  This  pulse  too  is  much  distorted  by  pickup. 

C.5  Conclusion 


In  conclusion,  it  may  be  said  that  the  exploding  wire  light 
source  has  been  developed  into  a  very  useful  and  reliable  piece  of  equipment. 
Used  in  conjunction  with  a  Kerr  Cell,  the  exploding  wire  source  has  been  the 
means  of  substantially  improving  the  time  resolution  of  the  interferograms 
taken  with  the  UTIAS  9"  Mach-Zehnder  interferometer.  Sample  interferograms 
taken  with  the  exploding  wire  source  are  shown  in  Figs.  2  and  3b,  and  the  im¬ 
proved  time  resolution  is  apparent  when  Figs.  2  and  3b  are  compared  with  Figs. 
1  and  3a. 


Upper:  5  volts/ div. 

1  psec. / div. 
Lower:  2  volts/div. 

0. 5  usee. /div. 


FIG. Cl  -  INTENSITY  PROFILE  OF  LIGHT  FROM  MAGNESIUM  SPA  TK 
LIGHT  SOURCE  -  TAKEN  WITH  FG&G  SD-100  PHOTODIO,  E 


Upper:  20  volts/div. 

5  usee. /div. 
Lower:  10  volt  /div. 
2  usee.  /div. 


FIG.C2  -  INTENSITY  PROFILE  OF  LIGHT  FROM  EXPLODING  WIRE 
LIGHT  SOURCF  -  TAKEN  WITH  EG&G  SD-100  PHOTODIODE 


FIG.  Cl  -  VIEWS  OF  EXPLODING  WIRE  LIGHT  SOURCE  INSTALLED 
IN  INTERFEROMETER 


FIG.C4  -  (  LOSE -TP  VIEWS  OF  I.ICIIT  SOI  KOI- 


PIG.C7  -  ELECTRICAL  CIRCUITRY  FOR  EXPLODING  WIRE  LIGHT  SOURCE  AND  KERR  CELL 
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a)  4680A  Hor:  Ipsec./div. 

Ver:  100  mv/div. 


b)  52loA  Hor;  Ip  c.  /div. 

Ver:  100  nv/div.- 


c)  5675A  Hor:  Ijtsec./div. 

Ver:  500  mv/div. 


FIG.  C ' 9  -  MONOCHROMATIC  INTENSITY  PROFILES  OF  LIGHT  FROM 

EXPLODING  WIRE  SOURCE,  TAKEN  AT  VARIOUS  WAVELENGTHS 
(Differences  in  output  reflect  variations  in  spectral  sensitivity  of 
the  photodiode  and  in  the  transmission  of  the  interference  filters) 


a)  Spectrogram  of  light  from  exploding  wire  source  viewed  directly, 
with  Kerr  Cell  removed.  Variations  in  intensity  in  the  spectrum 
reflect  variations  in  spectral  sensitivity  of  the  film  as  well  as 
variations  in  the  spectral  content  of  the  light. 
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b)  Spectrogram  of  light  from  exploding  wire  source  viewed  through 
the  Kerr  Cell  with  the  Kerr  Cell  cocked  open,  showing  the  cutoff 
due  to  the  Nitrobenzene  of  the  Kerr  Cell. 
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c)  Spectrogram  of  light  from  exploding  wire  source  viewed  through 
Kerr  Cell,  with  Kerr  Cell  operating,  i.e.,  open  for 
200  nanoseconds  only. 


FIG.  CIO  -  SPECTROGRAMS  OK  LIGHT  FROM  1  XPI.OD1NG  WIKI-  SOURCE 
UNDER  VARIOUS  CONDITIONS 

(Taken  with  a  Ililger  Medium  Quartz  Spectr  ograph  using 
Ilford  HI’S  Film.  Calibrating  spectrum  is  Mercury) 
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Upper:  l  psec. /dlv. 

10  volts/div. 
Lower:  l  psec.  /div. 

500  mv/div. 


FIG.Cll  -  TYPICAL  OSCILLOSCOPE  TRACE  USED  TO  CHECK  SYNCHRONIZATION 
BETWEEN  EXPLODING  WIRE  AND  KERR  CELL 
Upper  trace:  Sample  (  30  volts  approx. )  of  35kv  pulse  which  operates 
Kerr  Cell 

Lower  trace:  Output  from  photodiode  placed  inside  interferometer 

housing  to  monitor  light  output  from  exploding  wire  source 
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